(3.\  MBBSS 
HHMt 


Digitized  by  the  Internet  Archive 
in  2019  with  funding  from 
University  of  Alberta  Libraries 


h  tt  ps  ://arc  h  i  ve .  o  rg/d  eta  i  I  s/S  a  I  e  vs  ky  1 979 


V 


THE  UNIVERSITY  OF  ALBERTA 


NAME  OF  AUTHOR 
TITLE  OF  THESIS 


RELEASE  FORM 
Fred  C.  Salevsky 

Vortices  in  the  Ground  State  of  the 
Two-dimensional  Spin  k  XY  model 


DEGREE  FOR  WHICH  THESIS  WAS  PRESENTED 
YEAR  THIS  DEGREE  GRANTED  . *???... 


M.Sc. 


Permission  is  hereby  granted  to  THE  UNIVERSITY  OF 
ALBERTA  LIBRARY  to  reproduce  single  copies  of  this 
thesis  and  to  lend  or  sell  such  copies  for  private, 
scholarly  or  scientific  research  purposes  only. 

The  author  reserves  other  publication  rights,  and 


neither  the  thesis  nor  extensive  extracts  from  it  may 
be  printed  or  otherwise  reproduced  without  the  author's 
written  permission. 


THE  UNIVERSITY  OF  ALBERTA 


VORTICES  IN  THE  GROUND  STATE  OF  THE  TWO-DIMENSIONAL 

SPIN  h  XY  MODEL 

by 

FRED  C.  S ALEVS KY 


A  THESIS 

SUBMITTED  TO  THE  FACULTY  OF  GRADUATE  STUDIES  AND  RESEARCH 
IN  PARTIAL  FULFILMENT  OF  THE  REQUIREMENTS  FOR  THE  DEGREE 

OF  MASTER  OF  SCIENCE 

IN 

THEORETICAL  PHYSICS 

DEPARTMENT  OF  PHYSICS 


EDMONTON ,  ALBERTA 


FALL,  1979 


. 


. 


THE  UNIVERSITY  OF  ALBERTA 


FACULTY  OF  GRADUATE  STUDIES  AND  RESEARCH 


The  undersigned  certify  that  they  have  read,  and 
recommend  to  the  Faculty  of  Graduate  Studies  and  Research, 
for  acceptance,  a  thesis  entitled  VORTICES  IN  THE  GROUND 
STATE  OF  THE  TWO-DIMENSIONAL  SPIN  h  XY  MODEL  submitted  by 
FRED  C.  SALEVSKY  in  partial  fulfilment  of  the  requirements 
for  the  degree  of  Master  of  Science  in  Theoretical  Physics. 


It  is  not  at  all  true  that  the  scientist 
goes  after  the  truth.  It  goes  after  him. 


Soren  Kierkegaard 


ABSTRACT 


The  study  of  phase  transitions  in  two-dimensional 
systems  which  possess  more  than  one  degree  of  freedom  has  been 
propelled  by  apparently  contradictory  evidence  for  the  existence 
of  such  transitions.  The  situation  was  clarified  by  the  work 
of  Kosterlitz  and  Thouless  who  introduced  the  idea  of 
topological  long-range  order  to  replace  the  more  familiar 
concept.  The  phase  transition  in,  for  example*  the  classical 
(s=oo)  XY  or  plane  rotator  model  is  then  viewed  as  correspond¬ 
ing  to  the  unbinding  of  pairs  of  topological  excitations 
(vortices  and  antivortices) .  Such  a  picture  has  inspired  much 
theoretical  work  and  has  been  experimentally  confirmed.  This 
thesis  is  the  beginning  of  a  study  of  the  nature  of  the  phase 
transition  in  the  two-dimensional  spin  \  XY  model  to  see  if 
this  transition  can  also  be  described  in  the  Kosterlitz- 
Thouless  (KT)  picture.  Vortex  and  vortex-antivortex  pair 
operators  appropriate  for  this  model  on  both  the  square  and 
triangular  lattices  have  been  constructed.  The  expectation 
values  of  these  operators  in  the  ground  state  are  obtained 
using  the  finite  cell  method  of  Oitmaa  and  Betts.  This 
calculational  method  is  described  in  some  detail.  The 
conclusion  of  this  work  is  that  although  some  vortices  and 
antivortices  are  present  in  the  ground  state  of  the  two- 
dimensional  spin  h  XY  model,  it  appears  that  all  these 
excitations  are  paired  and  that  no  free  vortices  are  present 
in  this  state.  The  interpretation  of  the  results  on  the 
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triangular  lattice  is  complicated  by  the  relatively  high 
density  of  vortices  in  the  ground  state  on  this  lattice. 
These  results  are  in  agreement  with  the  KT  theory. 
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CHAPTER  I 


INTRODUCTION 

A.  Phase  Transitions  in  Two-Dimensional  Systems 

Convincing  arguments  exist  to  show  that  long-range  order 
is  not  present  at  non-zero  temperature  for  a  wide  variety  of  two- 
dimensional  systems.  As  early  as  1935,  Peierls  {1.1}  argued 
that  thermal  motion  of  long-wavelength  phonons  would  destroy 
the  long-range  order  of  a  two-dimensional  solid  in  the  sense 
that  the  mean  square  deviation  of  an  atom  from  its  equilibrium 
position  increases  logarithmically  with  the  size  of  the  system, 
and  the  Bragg  peaks  of  the  diffraction  pattern  formed  by  the 
system  are  broad  instead  of  sharp.  Mermin  showed  later  {1.2}, 
through  the  application  of  rigorous  operator  identities,  that 
long-range  order  of  this  simple  type  is  indeed  absent. 

Application  of  similar  arguments  can  be  made  to  show  that 
there  is  no  spontaneous  magnetization  in  a  two-dimensional 
magnet  with  spins  with  more  than  one  degree  of  freedom  {1.3}. 

This  latter  condition  allows  for  the  well-known  phase 
transition  in  the  two-dimensional  Ising  model.  It  is  also 
known  {1.4}  that  the  expectation  value  of  the  superfluid  order 
parameter  in  a  two  dimensional  Bose  fluid  is  zero. 

Contrasted  with  these  results  is  persuasive  evidence, 
from  various  sources,  for  some  kind  of  phase  transition  in 
two  dimensional  systems.  These  sources  include  direct  numerical 
work  on  a  system  of  hard  discs  {1.5}  and  high  temperature 
series  expansions  {1.6},  which  predict  a  diverging  susceptibility. 
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The  evidence  for  such  a  transition  is  much  stronger  for  the 
XY  model  than  for  the  Heisenberg  model {1.7}.  While  these 
results  apply  to  the  plane  rotator  model,  which  is  equivalent 
to  the  S=oo  XY  model,  a  phase  transition  is  also  believed  to 
occur  for  the  extreme  quantum  case,  the  spin  %  XY  model  {1.8}. 
If  long  range  order  as  conventionally  defined  is  not  present  at 
non-zero  temperature,  what  is  the  nature  of  the  predicted 
phase  transition? 

The  physical  picture  which  has  emerged,  largely  in  the 
works  of  Kosterlitz  and  Thouless  {1.9},  is  one  in  which 
"topological"  excitations,  in  the  form  of  vortices  and 
antivortices  for  the  plane  rotator,  play  the  primary  role. 

This  new  type  of  long-range  order,  of  a  topological  nature, 
may  exist  for  a  two-dimensional  solid,  neutral  superfluid 
and  for  the  XY  model,  but  not  for  a  superconductor  or  isotropic 
Heisenberg  model.  In  the  case  of  a  two  dimensional  solid, 
the  disappearance  of  topological  long-range  order  is  associated 
with  a  transition  from  an  elastic  to  a  fluid  response  to  a 
small  external  shear  stress. 

The  definition  of  long-range  order  which  is  adopted  in 
this  situation  arises  naturally  from  the  dislocation  theory 
of  melting  {1.10},  where  it  is  supposed  that  a  liquid  close 
to  its  freezing  point  has  a  local  structure  similar  to  that  of 
a  solid.  In  the  equilibrium  configurations  of  the  liquid, 
there  is  some  concentration  of  dislocations  which  can  move  to 
the  surface  under  the  influence  of  an  arbitrarily  small  shear 
stress  and  so  produce  viscous  flow,  while  in  the  solid  state 
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there  are  no  free  dislocations  in  equilibrium,  and  so  the 
system  is  rigid.  Now  although  isolated  dislocations  cannot 
occur  at  low  temperatures  in  a  large  system,  since  their 
energy  increases  logarithmically  with  the  size  of  the  system, 
pairs  of  dislocations  with  equal  and  opposite  Burgers  vector 
have  finite  energy  and  are  created  by  thermal  excitation. 

These  pairs  can  respond  to  an  applied  stress  and  so  reduce 
the  rigidity  modulus.  At  sufficiently  high  temperatures, 
the  largest  pairs  become  unstable  under  an  applied  shear  stress 
and  produce  a  viscous  response  to  the  shear.  The  presence 
or  absence  of  free  dislocations  can  be  identified  then  with 
the  absence  or  presence  of  topological  long-range  order.  It 
should  be  mentioned  at  this  point  that  the  theory  of  melting 
in  such  two-dimensional  crystals  and  related  problems  of  liquid- 
crystal  phases  in  very  thin  films  has  been  the  subject  of 
much  recent  activity  {1.11},  {1.12}. 

The  topological  excitation  in  the  plane  rotator  model  , 
or  neutral  superfluid,  which  is  analogous  to  the  dislocation 

in  a  two-dimensional  crystal  is  the  vortex.  The  plane 

rotator  Hamiltonian  is 

H  =  -J  E  S.-S.  =  -J  E  cos  (<b . -cj>  . )  ,  (1.1) 

<i,j>  1  3  <i,j>  1  3 

where  J>0,  II  S  Jl  =  1  and  <|k  is  the  angle  between  and  some 
arbitrary  reference  direction.  The  symbol  <i,j>  indicates 
that  the  sum  is  taken  once  over  every  bond  pair  (i,j).  In 
the  superfluid  case,  (j)  represents  the  phase  of  the  complex 
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superfluid  order  parameter.  By  expanding  the  argument  of  the 

2 

cosine  function  to  order  (<j>^-<j>j)  and  taking  a  continuum  limit 
it  is  possible  to  show  (1.9)  that  the  energy  of  an  isolated 
vortex  takes  the  form 


E  -  E  -  TrJln  (R/S )  ,  (1.2) 

where  R  is  the  size  of  the  lattice,  6  is  the  lattice  spacing 
and  Eq  arises  in  the  first  term  of  the  expansion  of  cos(<f>^“$j) 
and  represents  the  energy  of  the  zero  temperature  state  of  the 
plane  rotator  (see  also  {1.19}).  Since  isolated  vortices  (and 
anti vortices)  have  a  logarithmically  divergent  energy  in  the 
thermodynamic  limit,  R/6->«>,  they  cannot  occur  at  very  low 
temperatures.  Now  the  entropy,  S ,  associated  with  a  free 
vortex  is  approximately 


S  -  2k  In (R/ 5 ) 
v  B 


(1.3) 


which  is  arrived  at  by  noting  that  the  vortex  center  can 

2 

occupy  N-(R/6)  lattice  sites.  Combining  (1.2)  and  (1.3), 

the  free  energy  associated  with  an  isolated  vortex, 

F  =  (E  -E  )  -  TS  ,  can  be  written  as 
v  vo  v 

F  -  (ttJ  -  2kT)  ln(R/<5)  ,  (1.4) 

V  B 


and  from  this  it  can  be  seen  that  when  the  temperature  is 
higher  than 


T 

c 


TTJ 

2kB 


/ 


(1.5) 
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the  total  free  energy  can  be  lowered  by  creating  a  vortex. 

In  the  theory  of  Kosterlitz  and  Thouless,  the  temperature 

T  separates  two  distinct  phases  of  the  system.  The  "high" 
c 

temperature  phase,  T>Tc,  contains  isolated  vortices  and 

antivortices,  while  below  T  there  are  no  isolated  vortices 

c 

but  closely  bound  vortex-antivortex  pairs.  The  phase 

transition  at  T  corresponds  to  the  unbinding  of  these 

c 

pairs  and  the  creation  of  free  or  isolated  vortices. 

That  the  phase  of  the  system  at  temperatures  below  T 
actually  contains  vortex-antivortex  pairs  follows  from  the 
result  that  such  pairs  have  an  attractive  potential  energy 
which  cancels  the  logarithmic  divergence  (1.2)  of  a  single 
vortex.  This  can  be  understood  in  terms  of  an  analogy 
between  the  spin  system  containing  a  number  of  vortices  and 
a  two-dimensional  plasma,  originally  made  by  Kosterlitz  and 
Thouless  (1.9,  esp.  (1973)}.  If  vorticity  is  identified  with 
electric  charge,  vortices  interact  in  the  same  way  as  two- 
dimensional  electric  charges,  which,  in  turn,  interact  as 
parallel  lines  of  charge  in  three  dimensions.  The  interaction 
potential  then  depends  logarithmically  on  the  separation  of 
the  vortices.  Just  as,  in  two-dimensional  electrostatics, 
a  single  charge  has  infinite  electric  field  energy  and  a 
dipole  (a  pair  of  equal  and  opposite  charges)  has  finite  field 
energy,  assuming  the  charges  have  finite  size,  so  an  isolated 
vortex  has  infinite  energy,  while  a  vortex-antivortex  pair 
has  finite  energy.  The  vortex  "dipoles"  have  a  relatively 
short  range  effect  on  spin  correlations  compared  with  isolated 
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vortices  (and  antivortices)  for  the  same  reason  that  electric 
dipoles  have  short-range  electric  fields  compared  with  those  of 
isolated  electric  charges.  Because  of  this  the  spin  correla¬ 
tion  function  T  (r)  has  the  same  power  law  form,  for  ||  r ||/5 >>1 , 
as  predicted  from  spin-wave  theory  ignoring  the  effects  of 

vortices,  in  the  phase  below  T  {1.13}.  Once  isolated  vortices 

o 

can  appear,  they  have  long  range  disordering  effects  on  the 
spins  and  are  expected  to  lead  to  an  exponential  decay  of  the 
correlation  function  T  (r)  ,  for  large  Hr  l|/5.  The  transition 
from  one  type  of  behaviour  for  T (r)  to  the  other  indicates 
the  existence  of  a  phase  transition. 

All  the  above  arguments  can  be  made  more  quantitative, 
using  renormalization  group  methods  {1.14}.  In  particular, 
it  is  easy  to  see  that  the  estimate  (1.5)  for  the  critical 
temperature  is  an  upper  bound.  This  critical  temperature 
has  been  calculated  as  the  temperature  at  which  a  pair  of 
topological  excitations  will  dissociate,  ignoring  the  effect 
of  other  pairs  in  the  system.  These  other  pairs  will  relax 
in  the  field  of  the  first  pair,  thereby  renormalizing  downward 
the  e.g.  rigidity  modulus,  in  a  two-dimensional  crystal,  and 
consequently  reducing  the  critical  temperature.  This  tempera¬ 
ture  is  given  {1.14}  by  the  solution  of  the  relation 


7TJ 


2k  T 
B  c 


-  1  =  2irexp  (- 


,2j 


2k_,T 
B  c 


) 


(1.6) 


Other  predictions  can  be  made  of,  for  example,  the  critical 

exponent  n  (T) ,  of  the  spin  correlation  V (r)  for  T<T  and 

c 
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magnetic  susceptibility,  X/  in  the  magnetic  realization  of 

the  planar  model.  This  latter  quantity  is  predicted  to 

diverge  as  T  approaches  T  from  above  and  remains  infinite 

c 

for  all  T<T  . 

c 

The  significance  of  the  Kosterlitz-Thouless  picture 
is  that  it  shows  that  a  phase  transition  can  be  driven  by  the 
combining  or  breaking  up  of  pairs  of  topological  excitations 
without  the  appearance  of  conventional  long  range  order  like 
that  of  a  ferromagnet.  Similar  topological  arguments  have 
attracted  increasing  attention  in  other  areas  of  physics, 
notably  the  study  of  elementary  particles  where  they  may  play 
a  part  in  understanding  the  fundamental  constituents  of  matter. 
With  such  an  appealing  view  of  the  phase  transition  in  the 
plane  rotator  model,  it  is  natural  to  wonder  whether  these 
same  considerations  are  applicable  to  the  predicted  phase 
transition  in  the  spin  ^  XY  model. 

B .  Vorticity  Operators  and  the  Finite  Cell  Method 

The  aim  of  this  thesis  is  to  begin  the  investigation 
of  the  nature  of  the  phase  transition  in  the  spin  h  XY  model 
by  examining  the  vorticity  of  the  ground  state  of  this  model. 
The  quantum  nature  of  this  model  introduces  complications 
in  two  areas.  First,  a  notion  of  vorticity  which  is 
appropriate  for  quantum  operators  must  be  established. 

This  will  be  done  by  constructing  an  operator,  v,  among  the 
eigenstates  of  which  is  included  a  vortex  (and  an  antivortex) 
configuration.  This  operator  will  not  commute  with  the 
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Hamiltonian  of  the  model  so  the  ground  state  will  not  be  an 
eigenstate  of  this  vorticity  operator,  but  the  expectation 
value  of  this  operator  in  the  ground  state  will  provide  an 
indication  of  the  mean  number  of  vortices,  per  lattice  site, 
in  this  state. 

The  second  problem  arises  because  of  the  complex  nature 
of  the  ground  state  of  the  spin  %  XY  model.  This  state  is 
not  yet  known  exactly  for  the  two-dimensional  lattices  which 
will  be  considered  here,  so  it  is  necessary  to  use  indirect 
means  to  determine  the  required  ground  state  properties.  The 
finite  cell  method,  as  outlined  by  Oitmaa  and  Betts  {1.15}, 
will  be  employed  in  these  calculations.  The  basis  of  this 
method  is  to  calculate  exactly  the  ground  state  and 
ground  state  spin  correlations,  for  a  sequence  of  small 
(N£19)  finite  lattices  of  N  spins  and  then  extrapolate  the 
results  to  the  infinite  lattice  case.  Such  a  procedure  has 
been  applied  to  the  calculation  of  the  ground  state  energy 
and  transverse  and  longitudinal  magnetizations  of  both  the 
spin  %  XY  magnet  and  the  Heisenberg  antif erromagnet  on  square 
and  hexagonal  lattices  {1.15}.  More  recently  calculations 
have  been  completed  yielding  the  ground  state  energy,  as  a 
function  of  the  anisotropy  parameter,  for  the  anisotropic 
Heisenberg  antif erromagnet  on  the  square  lattice  {1.16}. 

Both  the  ground  state  energy  and  transverse  magnetization 
for  the  spin  h  XY  magnet  on  the  triangular  lattice  are  now 
also  known  from  such  finite  cell  calculations  {1.17}. 
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In  the  next  chapter  the  finite  cell  calculations  are 
described  in  some  detail.  The  original  results  of  Oitmaa  and 
Betts  (see  also  {1.18}),  for  the  ground  state  energy  of 
the  spin  H  XY  model  on  the  square  lattice, have  been  reproduced. 
While  this  result  was  obtained  using  a  sequence  of  cells 
with  an  even  number  of  spins,  further  calculations  of  this 
energy  for  several  cells  containing  an  odd  number  of  spins 
are  presented  as  well.  In  the  third  chapter,  the  concept  of 
vorticity  is  introduced  and  quickly  specialized  to  the  case 
of  a  vector  field  defined  on  a  lattice.  The  vortex  operators 
for  the  square  and  triangular  lattices  are  constructed,  as 
are  operators  which  are  relevant  to  the  study  of  vortex- 
antivortex  bound  pairs  in  the  ground  state  of  the  spin  h  XY 
model.  Once  these  operators  are  known,  their  expectation 
values  can  be  expressed  in  terms  of  spin  correlations 
available  from  Chapter  II.  These  expectation  values  are 
presented  and  discussed  in  the  last  section  of  Chapter  III, 

The  final  chapter  presents  a  summary  of  the  results  obtained 
and  further  discussion  on  the  implications  of  some  of  these 
results.  Suggestions  for  further  research  are  included  here 
and  the  section  ends  with  a  brief  discussion  of  experimental 
realizations  of  the  ideas  contained  in  this  thesis.  Three 
small  appendices  are  included  at  the  end  of  this  thesis.  The 
first  is  a  compilation  of  two  spin  correlations  for  finite 
cells  on  the  square  lattice,  while  the  second  is  a  brief 
description  of  the  contents  of  two  computer  files  which  were 
generated  during  the  course  of  this  work.  The  final  appendix 
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contains  a  construction  of  the  ground  state  of  the  anti¬ 
ferromagnetic  planar  model  on  the  triangular  lattice  which 
was  obtained  following  consideration  of  vorticity  on  the 
triangular  lattice.  As  mentioned  above,  this  thesis  is  but 
a  part  of  the  study  of  vortices  in  the  spin  h  XY  model.  High 
temperature  series  expansions  of  the  vortex  and  vortex- 
antivortex  pair  operators  derived  in  chapter  three  are  also 
being  carried  out.  The  complexity  of  such  calculations 
however  is  restricting  analysis  to  the  triangular  lattice 
at  present.  So  although  the  applicability  of  finite  cell 
calculations  to  this  problem  may  be  limited,  it  is  certain 
that  knowledge  of  the  nature  of  the  phase  transition  in  the 
two-dimensional  spin  h  XY  model  is  only  beginning  to  accumulate. 


CHAPTER  II 


GROUND  STATE  PROPERTIES  OF  THE  SPIN  h  XY  MODEL 

A.  Magnetization  of  the  Ground  State 

The  problem  to  be  examined  in  this  chapter  is  that 
of  determining  several  ground  state  properties  of  the  spin 
h  XY  model,  described  by  the  Hamiltonian 

=  -  2J  E  (S*S*  +  sYsY)  (2.1) 

XY  <i,j>  1  3  1  3 

X  V 

with  J>0.  The  operators  S^,  Si  are  components  of  a  spin  h 
operator  S^  associated  with  each  lattice  site  i  and  the 
symbol  <i,j>  indicates  that  the  sum  is  taken  over  nearest 
neighbour  pairs  of  sites  on  a  regular  two-dimensional  lattice. 
The  ground  state  energy  Eq  e  <0|Hxy|0>  will  be  the  focus  of 
this  chapter  but  it  is  obvious  that  Eq  is  proportional  to  the 

X  X 

nearest  neighbour  correlation,  E  =  -  2qNJ<0 | SQS ^ | 0> ,  where 
6  is  the  lattice  spacing  and  q  the  coordination  number  of  the 
lattice.  The  calculational  method  can  easily  be  extended  to 
yield  correlations  between  nth-nearest  neighbour  spin 
components . 

Unlike,  for  example,  the  Ising  model 

H  =  -  J  £  SZSZ  ,  (2.2) 

1  .  i  1 

<i , 3  >  J 

for  which  the  ground  state  is  clearly  that  in  which  the  spins 
on  all  lattice  sites  are  parallel,  the  ground  state  of  H XY 
cannot  be  written  down  by  inspection  and  in  fact  is  not 
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known  exactly,  in  two  or  three  dimensions.  Rather  than 
attempting  an  exact  solution  of  the  problem  on  an  infinite 
lattice,  the  basic  idea  of  the  following  calculations  has 
been  to  carry  out  exact  determinations  of  ground  state  energy 
(or  correlations)  for  a  sequence  of  small  finite  lattices  of 
N  spins  and  then  to  extropolate  to  the  infinite  lattice  case. 
The  N-spin  cells  are  endowed  with  periodic  boundary  conditions 
so  that  an  infinite  lattice  can  be  regarded  as  tiled  with 
N-spin  units.  The  limit  N->°°  (or  N  ^-*-0)  then  produces  results 
appropriate  for  the  infinite  lattice  "tiled"  with  a  single 
unit. 

For  a  finite  spin  h  system  it  is  in  principle  a  straight 
forward  matter  to  find  the  energy  eigenvalues,  the  problem 

N 

reducing  to  the  diagonalization  of  the  Hamiltonian  in  the  2 
dimensional  space  of  states.  Such  a  procedure  becomes 
computationally  intractable  around  N=8  unless  symmetry 
considerations  can  be  used  to  block  diagonalize  the 
Hamiltonian  matrix.  To  this  end  consider  the  longitudinal 
magnetization  operator 

M  =  Z  SZ  ,  (2.3) 

z  i  ! 

which  commutes  with  H  and  thus  provides  a  quantum  number  to 

A  x 

label  the  eigenstates  of  the  Hamiltonian.  The  following 

proof,  due  to  Mattis  (2.1),  shows  that  in  the  ground  state 

$  ,  the  quantum  number  m  is  zero.  It  is  first  convenient 
o  z 

to  perform  the  rotation  in  spin  space  which  takes  to 

and  S?  to  -sY,  at  all  lattice  sites  i.  Then  the  operators 
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H  and  M  have,  in  this  new  representation,  the  form 

Ax  Z 


H 


XY 


2  J  Z 
<i,  j> 


(SXSX  +  szsz) 

13  1 


(2.4) 


and  M  = 

z 


Z  S*  =  h  i  Z  (s;  -  s. ) 
i  1  i  1  1 


(2.5) 


where  the  raising  and  lowering  operators  have  been 
introduced  according  to  the  definitions 


II 

+  • 
VI 

sx  + 

isT 

1 

1 

1 

sT  = 

1 

sx  - 
1 

isY 

1 

• 

(2.6) 

Since  the  rotation  in  spin  space  is  a  similarity  transforma¬ 
tion,  the  eigenvalues  of  H  and  M  remain  the  same;  in 

AX  Z 

particular,  the  value  of  m  when  the  energy  is  minimum  will 

z 

be  unchanged  by  this  transformation.  The  Hilbert  space  of 
states  for  a  lattice  with  N  sites  consists  of  2N  configurations 
of  two  distinct  types: 


+  Pn  , 

l/j  =  C  7T  (S  )  0> 

ra,ev.  a  n  1 

n 

+  Pn, 

and  ip  ,  =  c  tt(S)  0> 
ra,od.  a  n  1 

n 


with  Zp  =  0 , 2 , 4 ,  .  .  .  ,  ( 2 . 7a) 
n 

with  Zp  =  1 , 3 , 5 , . . . , ( 2 . 7b) 
n 


where  | 0>  is  the  state  with  all  spins  "down",  c^  is  a 
normalization  constant  and  the  set  of  integers  {p^}  defines 
the  state  By  using  (2.6)  for  in  (2.4)  it  can  be 

seen  that  has  no  matrix  elements  to  connect  the  "even" 

and  "odd"  states.  The  two  subspaces  are  decoupled  and  the 
ground  state  of  each  must  be  studied  individually. 


r 
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2 

In  the  representation  (2.4)  and  (2.5)  the  are  all 

diagonal,  but  the  operators  are  not.  The  ground  state  in 

either  even  or  odd  subspace  takes  the  form 

*  „  =  E  F ( r)  ^  ,  with  Z|F(r)  |2  =  1  ,  (2.8) 

v o,r  a  Ka,r  1  a  1 

a  a 

and,  according  to  a  well-known  theorem  of  Frobenius  (2.2), 

( a) 

has  the  property  that  all  the  F  ,  for  r  =  even  or  odd, 

( IT ) 

can  be  chosen  real  and  positive.  None  of  the  Fa  vanishes, 
since  (2.4)  connects  all  configurations  within  either  sub- 

(  27  ) 

space,  nor  is  any  F  of  opposite  (negative)  sign  in  the 

06 

ground  state.  The  proof  of  this  statement  is  by  contradiction. 
If  some  amplitude  were  not  positive,  the  variational  energy 
e  °  =  <$Q  r  |  |  r>  could  be  decreased  by  making  it  so. 

However,  e°  is  already  the  lowest  possible  energy  for  the 

(  2T  ) 

respective  subspace,  and  hence  all  F  >0.  Finally,  the 

06 

ground  state  in  each  subspace  is  non-degenerate,  as  no  other 

eigenstate  of  H v  can  satisfy  the  condition  of  having  all 

positive  amplitudes  and  yet  be  orthogonal  to  r* 

Consider  first  the  state  §  .  Since  M_  as  given 

o,even  Z 

by  (2.5)  connects  the  even  and  odd  subspaces,  and  in  fact 

takes  even  states  into  odd  states,  the  only  way  $  can 

1  1  o,ev. 

be  a  simultaneous  eigenstate  of  and  Mz  is  if 


M„  4 
Z  o,ev 


0 


(2.9) 


This  argument  also  applies  to  so  that  the  ground  state 

magnetization  of  the  XY  magnet  is  zero.  Thus,  it  suffices  to 

N 

examine  the  niz=0  subspace  of  the  2  dimensional  Hilbert  space 


of  states. 


-  I 
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Actually  a  complication  arises  if  N,  the  number  of 

sites  per  cell,  is  odd.  In  this  case  there  is  not  state  in 

either  the  even  or  odd  subspace  which  has  m  =0,  the  minimum 

z 

value  of  |m  |  being  By  rotating  (2.4)  and  (2.5)  about  the 

z 

X  V  V  z  z 

S.  axis  in  spin  space  such  that  S.  -*-s  •  and  S .  -*-S .  ,  states 

l  i  i  l  l 

in  the  even  and  odd  subspaces  are  interchanged,  M  is  mapped 

z 

into  -M  and  H  is  left  invariant.  This  serves  to  illustrate 

Z  A1 

an  essential  degeneracy  of  the  two  subspaces  (even  and  odd) 
which  is  missing  when  N  is  an  even  integer.  Because  of  this 
degeneracy,  and  unlike  the  even  N  situation,  it  is  possible 
to  have  an  eigenstate  of  H  which  is  some  linear  combination 

X  JL 

of  eigenstates  from  both  even  and  odd  subspaces  and  which  is 

also  an  eigenstate  of  M  ,  with  non-zero  eigenvalue. 

z 

The  value  of  m  in  the  odd  N  ground  state  is  ±h, 

z 

depending  on  the  chosen  linear  combination  of  $  and 

r  o , ev . 

$o  To  prove  this  rigorously  consider  a  reference 

Hamiltonian 

Href*  =  -  2 J /N [ ( ISX ) 2  +  (IS2)2]  (2.10) 

n  n  J 

n  n 

in  which  every  spin  on  the  same  lattice  as  before  interacts 

with  every  other  spin.  This  is  a  molecular  field  limit  of 

the  Hamiltonian  (2.4).  On  the  one  hand  the  ground  states 

re  f 

in  the  even  and  odd  subspaces  (2.7),  r *  all  have  positive 
amplitudes,  by  Frobenius'  theorem.  They  are  therefore  not 
orthogonal  to  their  counterparts  in  Equation  (2.8),  and 
therefore  share  the  same  quantum  number  m  .  But  now  the 
energy  levels  of  (2.10)  are  easily  obtained  (by  back-rotating 
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s'r+sj')  as  e  =  -2J/N[l  (1+1)  -r/]  with  I=I.T,  I..-1,  1-2,..., 

n  n  ret.  z  N  N  N 

m  =1,  1-1 , . . . , -1+1 , -I  and  I  =N/2.  Evidently  the  ground  states 

^  IN 

of  H  .  belong  to  1=1  and  m  =0,  when  N  is  even,  or  m  =±% 
it  g  r  t  z  z 

when  N  is  odd. 

To  summarize  the  above  arguments,  it  is  sufficient  to 

note  that  when  considering  the  XY  model  on  a  lattice  of  N 

spin  V s ,  one  need  only  consider  the  m  =0  subspace  of  the  2N 

z 

dimensional  Hilbert  space,  when  N  is  an  even  integer.  When  N 

is  odd  it  is  sufficient  to  consider  either  the  m  =%  or  m  =-h 

z  z 

subspaces.  It  should  also  be  mentioned  that  these  results 

apply  in  more  generality  than  indicated  by  the  above  proof. 

For  example,  it  can  be  seen  that  nothing  requires  all  the  spins 

to  be  equal  and  in  fact  the  results  above  are  easily  generaliza- 

ble  to  the  case  when  the  total  spin  on  each  site  i  is  S^, 

the  place  of  N  is  taken  by  2ES..  Further  details  and 

i 

extensions  are  available  in  the  original  reference  {2.1}. 

B .  Cell  Geometry  and  Symmetry  Group 

The  above  effort  has  succeeded  in  reducing  the  number 

of  states  which  must  be  considered  in  diagonalizing  H XY  from 

2N  to  N! /[ (N/2) ! ]*  2 ,  for  N  even  or  N!/[ (N+l/2) ! (N-l) /2) ! ] ,  for 

N  odd.  These  latter  are  the  numbers  of  states  in  the  m  =0 

z 

and  m  =h  subspaces  respectively.  For  example,  there  are 

4i 

NI/(N/2)I  ways  of  placing  N/2  objects  ("up"  spins)  on  N 
lattice  sites;  the  extra  factor  (N/2) I  corrects  for  the 
indistinguishability  of  the  spins  placed  on  the  lattice. 

These  numbers  are  still  quite  large.  For  N=16  there  are 


17 


1 6 

12,870  states  in  the  m  =0  subspace  or  ^  20%  of  the  total  2 
states.  A  further  reduction  in  the  number  of  states  which 
must  be  considered  to  diagonalize  H  is  needed  to  be  able 
to  push  the  calculation  to  cells  with  N>10.  Such  a  reduction 
is  afforded  by  consideration  of  the  space  group  symmetry  of 
the  particular  finite  cell,  since  the  ground  state  eigenvec¬ 
tor  must  transform  according  to  one  of  the  irreducible 
representations  of  the  space  group.  The  same  Frobenius 
theorem  mentioned  above  requires  the  ground  state  eigenvector 
to  transform  according  to  the  identity  representation.  That 
is,  the  ground  state  is  invariant  under  the  action  of  each 
element  of  the  space  group. 

The  determination  of  the  space  group  relevant  for 
each  finite  cell  begins  with  the  specification  of  the  geometry 
of  the  cells  upon  which  the  calculation  will  be  based.  For 
the  square  lattice  cells  were  chosen  from  two  series:  an 
even  series,  with  cells  of  4,  8,  10,  16  and  18  spins  each, 
and  an  odd  series,  with  cells  of  4 ,  9  and  13  spins.  These  are 
shown  in  Figure  2.1.  Although  there  are  many  possible  choices 
for  both  the  number  of  spins  per  cell  and  the  geometry  of  a 
cell  with  a  given  number  of  spins,  the  choices  made  for  this 
calculation  are  such  that  the  symmetry  group  of  each  cell  is 
as  large  as  possible.  All  the  cells  considered  here  have  a 
square  border,  a  characteristic  which  increases  the  dimension 
of  the  point  group  (i.e.  rotations  and  reflections)  relative 
to  that  of  cells  with  rectangular  or  rhomboidal  boundaries. 

The  number  of  spins  which  can  be  included  in  such  square  cells 
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Figure  2.1  -  Cell  geometries  on  the  square  lattice. 
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2  2 

is  given  by  n  +m  ,  with  n,m  integer.  The  numbers  n  and  m 

represent  the  coordinate  distance  between  cell  centers  (i.e. 

/*2  ? 

3  cells  across,  1  cell  up)  and  so  /n  +m  represents  the 

distance  between  cell  centers.  This  distance  is  also  the 

2  2 

length  of  an  edge  of  the  square  boundary  and  so  n  +m  is 

2 

both  the  area  in  units  of  6  ,  where  5  is  lattice  spacing,  of 

a  cell  and  the  number  of  spins,  or  lattice  sites,  the  cell 

contains.  This  relation  predicts  an  additional  cell  belonging 

2  2 

to  the  odd  series,  with  N=17=(4)  +(1)  ,  but  calculations  were 

not  made  for  this  value  of  N. 

Once  the  cell  geometry  is  chosen,  an  infinite  lattice 
is  tiled  with  N-spin  cells  and  each  cell  is  numbered,  as 
shown  for  the  9-spin  cell  in  Figure  2.2.  Two  groups  are 
needed  in  order  to  be  able  to  generate  the  full  space  group 
S„.  The  first  is  the  group  of  translations,  Txt,  which  consist 
of  those  transformations  on  the  tiled  lattice  which  move  each 
of  the  N  spin  sites  to  a  given  fixed  site.  These  transforma¬ 
tions  can  be  visualized  as  moving  the  grid  of  cell  boundaries 
over  the  underlying  numbered  lattice.  The  second  is  PN,  the 

group  of  point  symmetries,  consisting  of  rotations  by  90°, 

2  3 

denoted  by  C^,  (C^)  and  (C^)  ,  and  reflections  about 

vertical,  horizontal  and  diagonal  axis  of  the  square  cell. 

Some  cells  do  not  appear  to  be  rotation  symmetric,  for  example 
the  N=8  and  10  cells.  However,  all  cells  do  have  this 
symmetry,  although  it  may  require  a  small  shift  of  the 
boundary  grid  over  the  lattice  to  reproduce  a  similar  pattern 
of  spins  in  each  cell.  What  is  exceptional  is  that  the  N=10 
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and  N=13  cells  are  not  reflection  symmetric,  the  point 
groups  S1Q  and  S^3  consisting  of  only  the  proper  rotations, 
including  the  identity. 

Using  the  groups  TN  and  PN  as  input,  a  computer  program 
(GROUNDS YMM) ,  available  in  the  Statistical  Physics  Program 
and  Data  Library)  was  used  to  generate  the  group  elements 
of  SN-  The  elements  of  these  transformation  groups  are 
stored  and  manipulated  as  permutations  of  N  numbers.  For 
example,  the  identity  element  of  S^  is  123456789, 
which  corresponds  to  sites  on  the  9- spin  cell  through  the 
numbering  illustrated  in  Figure  2.2.  Four  sample  symmetry 
transformations  are  illustrated  in  Figure  2.3.  Rotation 
about  the  center  of  symmetry  by  270  ,  (C^)  ,  is  represented 

as  a  permutation  by  369258147,  which  means  that 
under  this  transformation  site  3  moves  to  site  1,  site  6 
moves  to  site  2,  ...,  etc.  Also  illustrated  are  the  reflec¬ 
tion  about  the  1-9  axis,  ,  and  a  sample  translation,  t^. 
With  the  symmetry  transformations  represented  by  permutations, 
the  rules  for  multiplying  permutations  allow  the  quick 

3 

calculation  of  such  compound  transformations  as  (C^)  «t3: 


(e 


Note  that  (C^) 


3 


123456789) 

312645978 

369258147 

258147369 

3 

t3  ^  t3»(C4)  ,  and  in  general  p^t^  ^ 


f 
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Figure  2.2  -  Numbering  of  the  9-spin  cell. 
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Figure  2.3  -  Sample  symmetry  operations  on  the  9-spin  cell. 
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p.eP., ,  t.eT...  Actually  GROUNDSYMM  does  not  simply  calculate 
i  N  3  N 

the  product  set  {p^'t^}.  The  program  multiplies  each  element 
held  in  memory  by  every  other  such  element  and  quits  only 
when  this  procedure  fails  to  generate  a  new  symmetry  trans¬ 
formation.  The  number  of  elements  in  each  group  S  ,  for  N>5, 
is  shown  in  the  second  column  of  Table  2.1.  It  is  interesting 
to  note  that  in  all  cases  the  order  of  S  ,  0(N) ,  is  simply 
the  product  of  the  orders  of  T._  and  PXT:  there  are  N  trans- 
lations,  including  the  identity,  for  each  N-spin  cell  while 
there  are  eight  members  of  PN  (three  rotations,  four  reflec¬ 
tions  and  the  identity)  for  N^IO  or  13.  In  these  special 
cases  P^  consists  of  just  proper  rotations  (including  the  identity) 
thus  has  order  4.  This  result  is  true  since  the  translation 

group  Txt  is  a  normal  subgroup  of  SXT  and  so  SXT  =  {p.*t.}  = 

3  N  N  N  *1  j 

{tfa’P^}  and  the  order  of  SN  is  the  product  of  the  orders  of 

T_.  and  P__ .  {2.3}. 

N  N 

C .  Basis  States 

The  ground  state  which  is  to  be  constructed  from  states 

in  the  m  =0  subspace  of  the  Hilbert  space  (2.7),  symbolized 
z 

as  i  must  be  invariant  under  the  action  of  the  space 

group  S  .  If  this  ground  state,  for  each  N-spin  cell,  is 

(N) 

denoted  by  ,  the  required  invariance  may  be  expressed  as 

g.$^N)  =  $  ,  for  all  g.  e  S.T  .  (2.11) 

^i  o  o  l  N 

If  the  members  of  the  mz=0  subspace  of  the  Hilbert  space  of 

(N)  .  (N) 

states  are  denoted  by  y  ,  then  it  is  clear  that  if  y 

J  Aa  a 
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Table  2.1  -  Ground  state  energy  data  for  the  spin  h  XY  model 

on  the  square  lattice. 


#of  spins/cell 

order  of 

N 

#of  basis 

-E  /J 
o 

-E  /NJ 
o 

N 

0  (N) 

states 

4 

— 

2 

5.6569 

1.4142 

8 

64 

6 

9.3808 

1.1726 

10 

40 

9 

11.5216 

1.1522 

16 

128 

153 

17.9984 

1.1249 

18 

144 

398 

20.1439 

1.1191 

5 

— 

1 

6.0000 

1.2000 

9 

72 

5 

10.3470 

1.1497 

13 

52 

38 

14.6330 

1.1256 
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enters  into  $  ^  with  coefficient  a^  so  must  each  state 

o  a 

g-xi^  for  all  g.eS  .  But  each  g-X^N^  is  just  one,  though 

1  06  1  IN  1  06 


possibly  the  same,  state  x^  ^  in  H 


Only  those  states 


4^N^  in  H  ^  which  are  not  related  to  each  other  by  a 

symmetry  transformation  will  enter  into  with  different 

(N) 

coefficients  a  .  In  this  way  such  states  are  said  to  form 


cx 


(N) 


an  unsymmetrized  basis  for  4>q  .  The  symmetrized  basis  for 

,(N)  .  ,  (N)  ,  /  (N)  v  ,  (N)  ^  (N) 

$  is  40  =  1/n  Z  g.i p'  .  The  integer  n  is 

o  a,sym.  a  .  i  a  a 

1  .  (N) 

the  number  of  group  elements  g^  which  leave  ip  invariant. 

Since  there  are  0(N)/n^N^  distinct  states  X^N^  generated  by 
(N) 

4>  and  the  sum  above  is  taken  over  all  0(N)  group  elements, 
06 

the  factor  l/n^N^  is  necessary  to  reproduce  each  X^N^ 
exactly  once. 

To  generate  the  basis  states,  use  was  made  of  the 
program  GROUNDBASIS,  also  in  the  Statistical  Physics  Program 
and  Data  Library,  modified  of  necessity  for  each  size  of  cell 


Taking  the  symmetry  group  S^  as  well  as  one  initial  state 
(N) 

4>^  and  its  symmetry  number  nQ  as  input,  this  program 


generates  the  complete  set  of  basis  states  ip 


(N) 


a 


and  their 


symmetry  numbers  n  .  This  is  accomplished  by  first  over¬ 
turning  nearest-neighbour  pairs,  with  opposite  spin,  in  order 

(N) 

to  generate  a  new  state  ip £  which  is  not  related  by  symmetry 

.  (N) 

transformation  to  ip^  •  More  than  one  state  inequivalent  to 

4^N^  may  be  generated  by  flipping  just  one  pair  of  opposite 

spins;  if  so  this  state  is  also  stored,  along  with  its 

(N) 

symmetry  number.  The  program  then  moves  on  to  4*^  and  tries 
to  generate  new  basis  states  by  flipping  nearest  neighbour 
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pairs  of  opposite  spin  from  this  state.  Eventually  the  point 

will  be  reached  at  which  no  new  basis  states  are  generated 

as  these  spin  pairs  are  overturned.  The  program  will  keep 

trying  to  generate  new  states  however,  until  a  maximum  number 

(N) 

of  spin  pairs,  relative  to  i|)  ,  have  been  flipped.  This 

number  must  be  fed  into  the  program  through  the  input  file. 
Experience  has  shown  that  all  basis  states  will  be  generated 
if  this  number  is  >N/2.  For  example,  the  last  basis  state  of 

*N/ 

the  16  spin  cell  is  generated  after  6  pairs  of  spins  have 
been  flipped  relative  to  the  initial  state.  Execution  stops 
once  the  maximum  number  of  spin  pairs  have  been  overturned 
and  the  list  of  basis  states,  with  their  symmetry  numbers  and 
data  concerning  the  numbers  of  spin  pairs  flipped  is  written 
into  the  output  file. 

The  unsymmetrized  basis  elements  for  the  9-spin  cell, 
along  with  their  symmetry  numbers , are  presented  in  Figure  2.4. 
Actually,  as  mentioned  above,  these  states  are  just  the 
generators  of  the  basis  states,  which  can  be  written  formally 
as  the  symmetrized  sum 


(N) 

a , sym . 


0(N) 

E 

i=l 


,  (N) 

g .  ip 

a 


(2.12) 


a 

In  addition,  we  require  the  basis  states  to  be  normalized. 

Since  there  are  0(N)/n^N^  distinct  states  in  each  symmetrized 

state  ,  and  because  all  the  individual  states  in  H ^ 

ra , sym.  o 

are  normalized,  the  proper  normalization  factor  is  (n^N^/°(N)) 

vA 

and  the  symmetric,  normalized  basis  states  can  be  written  as 


h 
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Figure  2.4  -  Basis  states  of  the  9-spin  cell  and  their 

symmetry  numbers.  (  =spin  "up",  =spin  "down") 


t(N) 

a , sym. , norm. 


1 


=  y(N) 

a 


(n^N)  0(N) )h 


0(N) 

S 

i=l 


(N) 

a 


(2.13) 

(N) 

The  number  of  symmetrized  and  normalized  basis  states  ¥ 
for  each  N  is  shown  in  the  third  column  of  Table  2.1.  Notice 
that  this  number  grows  rapidly  with  increasing  N;  this  is  to 
be  expected  since  the  number  of  states  in  is  ( N I )  /  [  (N/2 )  !  ] 

for  even  N,  while  the  number  of  elements  in  SN  is  proportional 
to  N.  As  the  number  of  basis  states  grows  r  so  does  the 
computing  time  necessary  to  generate  them,  and  the  limit  of 
computational  feasibility  seems  to  reside  at  N=18  on  the  square 
lattice  and  N=19  on  the  triangular  lattice.  This  latter  cell 
has  838  basis  states. 


D.  Matrix  Elements  of  the  Hamiltonian  and  Ground  State  Energy 
Using  the  raising  and  lowering  operators  defined  in 
(2.6),  it  is  straightforward  to  recast  the  XY  model  Hamiltonian 
of  (2.1)  in  the  form 

hxy  =  -J  <.E.>  (sisj  +  sjsi)  (2-14) 

For  a  finite  cell  with  periodic  boundary  conditions,  every 
lattice  site  on  the  square  lattice  has  four  nearest  neigh¬ 
bours.  As  an  example,  see  Figure  2.2  where  the  dark  lines 
are  the  structural  links  of  the  original  cell  and  the  dotted 
lines  represent  links  due  to  the  periodic  boundary  conditions. 
These  can  be  thought  of  as  "hooks"  by  means  of  which  an 
infinite  lattice  is  constructed  with  tiles  of  9  spins  each. 
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In  the  Hamiltonian  (2.14)  each  nearest  neighbour  link  is 
to  be  counted  only  once.  The  number  of  such  distinct  nearest 
neighbour  links  is  2N. 

Thus  for  the  9  spin  cell  the  Hamiltonian  has  eighteen 
terms  and  can  be  written  in  the  following  shorthand  form: 


"HXY/J  =  (12)+(23)+(14)+(25)+(36)+(45)+(56)+(47)+(58) 

+  (69)+(78)+(89)+(17)+(28)+(39)+(13)+(46)+(79), 


(2.15) 


where  (ij)  =  S+S .  +  S . S+  . 

1  3  13 


It  is  the  matrix  elements  of  the  operator  -H^y/J, 
rather  than  H  ,  w^ich  are  calculated  by  the  program 
GROUNDXYHAM,  the  third  and  final  member  of  the  sequence 
available  in  the  Statistical  Physics  Program  and  Data  Library 
This  program  also  uses  a  standard  diagonalization  routine  to 
yield  the  eigenvalues  and  eigenvectors  of  ~Hyy/J  for  all 
N£l6.  For  N=18  the  iterative  power  method  {2.4}  must  be 
used  to  generate  the  largest  eigenvalue  and  corresponding 
eigenvector.  This  will  be  described  subsequently. 

For  the  moment  consider  an  operator  A ,  slightly  more 
general  than  the  Hamiltonian,  in  that  it  may  consist  of 
products  of  any  number  of  spin  operators  on  any  combination 


of  sites.  The  problem  of  evaluating  the  matrix  elements  of 


A  in  the  basis  ¥ 


(N) 


a 


is  that  of  evaluating  expressions  of  the 


form 
(N) 


3 


A  Y 


(N)>  = 


a 


1  „  ,  (N)  1  , 

(n)  (n).%  gj 
na  ne  >  J 


1  (Z  Ag.  |i()  *N)  >) 

0 (N)  • 


(2.16) 


I 
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Now  the  sum  Z  g.<^„ 

j  3  3 


(N) 


can  be  written  as  (£  gk<^N^  1)9-^ 

k 


since  in  both  cases  the  summation  is  over  all  elements  geS„. 

N 

Then  (2.16)  becomes 


<V(N)  |a|¥(N)>  = 
6  1  1  a 


,  (N)  (N),-j 

na  nB 


^  {l  V*r1|)(k  l  gi1*gi)  l+aN)> 

(2.17) 


Specializing  to  the  case  of  which  is  invariant  under 


the  action  of  S  ,  equation  (2.17)  reduces  to 


(N)  l  HXYk,  (N) 


3 


a 


>~ 


_ 1 _  iZ  q  <^(N)|) 

,  (N)  (N)  .3*  v  gk  ^3  1  ’  J  ' 

(n  nD  )  k 


a 


a 


•3 


(2.18) 


Operationally,  equation  (2.18)  has  the  following 

meaning.  The  action  of  H^y  on  the  basis  generator  is 

known  from  the  action  of  the  operators  (ij);  if  the  spins  on 

(N) 

sites  i  and  j  of  \p  are  different  they  are  flipped  and  if 
these  spins  are  the  same  the  state  is  annihilated.  What 
GROUNDXYHAM  does  is  examine  each  element  of  the  set  {g.  } 

K  p 
(N) 

to  find  all  the  states  which  differ  from  ip  at  two,  and  only 

Ya  J 

two,  lattice  sites.  This  ensures  that  the  two  spins  on  these 

sites  are  different  (one  "up"  and  one  "down").  The  program 

subsequently  checks  the  location  of  these  two  spins  and  counts 

each  of  those  states  for  which  the  two  differing  spins  are 

nearest  neighbours.  In  this  way  the  matrix  element 

(£  g,  I) — is  generated  for  given  a,  3.  This 

^  k  p  j  a 

element  is  always  an  integer.  To  further  reduce  computing 
time  only  those  ,  for  a  given  a /  need  be  considered  as  the 
operator  -H VV/J  is  an  observable  and  must  have  symmetric  matrix 

A  JL 
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elements.  The  final  step  in  the  program  is  the  multiplication 

of  the  raw  matrix  elements  by  the  normalization  factor 
,  (N)  (n)  -h 

(n^  n^  ' )  .  The  total  set  of  matrix  elements  (2.18)  is 

then  stored  awaiting  the  action  of  the  diagonalization 
routine . 

As  mentioned  above,  a  standard  routine  was  used  to 
calculate  the  complete  set  of  eigenvalues  and  corresponding 
eigenvectors  of  the  operator  -H^/J  for  8£N£16.  The  calcula¬ 
tions  for  N=4  and  5  were  completed  "by  hand".  The  resulting 
ground  state  energy  -E  /J  was  obtained  from  the  list  of 
eigenvalues  by  inspection,  and  both  -E  /J  and  -Eq/NJ  are 
tabulated  in  the  last  two  columns  of  Table  2.1.  For  N=18, 
the  matrix  of  -H  Y/J  is  398  by  398,  which  is  above  the  limit 
of  efficiency  and  accuracy  of  the  standard  diagonalization 
routines.  Since  only  the  largest  eigenvalue  and  corresponding 
eigenvector  are  required,  the  iterative  power  method  was 
used  to  generate  these  quantities.  An  initial  normalized 
vector  <J>  is  chosen  and  the  vector  <b  ’  =  J)<f>  is 

calculated.  This  vector  is  in  turn  normalized,  (J)  ^=<|>^/||<j)^  1 1 
and  the  procedure  repeated  with  <f>^.  After  n-25  iterations, 
the  vector  <j)  has  converged  to  the  ground  state  eigenvector 
while  ||  4)n'||/||4)nll  has  converged  to  the  maximum  eigenvalue  of 
-H  „r/J,  i.e.  -E  /J.  A  suitable  choice  for  the  initial  vector 

d>  is  one  in  which  all  entries  are  the  same,  i.e.  1,  with  an 
r  o 

-k  • 

overall  normalization  factor  of  e.g.  (398)  in  the  case  of 
the  18-spin  cell  Hamiltonian  matrix.  In  this  manner  the 
eigenvalue  -Eq/J=20 . 1439  was  obtained. 


I 


f 
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Figure  2.5  is  a  graph  of  the  ground  state  energy/spin, 
-Eq/NJ  vs.  1/N.  The  results  for  N  an  even  integer  lie  on  a 
straight  line,  while  the  odd  N  data  do  not  lie  on  this  line. 
Note  however  that  a  line  drawn  through  the  estimates  for  N=9 
and  13  is  parallel  to  the  former  line.  The  series  corresponding 
to  even  N  confirms  the  results  of  Oitmaa  and  Betts  {2.5}  and 
leads  to  a  value  of  ground  state  energy/site,  Eo/NJ,  in  the 
infinite  lattice  of  -1 . 078±0 . 002 .  Although  only  two  members 
of  the  odd  N  series  are  presented,  since  the  results  for  N=5 
do  not  lie  on  the  line  formed  by  the  larger  N  data  and  in 
fact  are  not  expected  to,  it  would  appear  that  the  value  of 
Eq/NJ  predicted  for  the  infinite  lattice  from  the  odd  series, 
-1.07±0.01,  is  not  consistent  with  the  even  N  extrapolation. 

The  larger  uncertainty  associated  the  odd  N  extrapolation 
reflects  the  use  of  only  two  data  points  and  certainly 
before  any  definite  conclusion  can  be  reached,  the  ground 
state  energy  for  the  N=17  cell  must  be  calculated.  More  will 
be  said  on  this  matter  in  Chapter  IV  -  Summary  and  Discussion. 


E .  Correlations 

Other  quantities,  besides  the  ground  state  energy, 
can  be  evaluated  once  the  basis  states  have  been  computed. 
Principal  among  these  quantities  are  two  and  four  spin 
correlations,  for  example 


<$<N> 

|sxsx|$(N)>  , 

o 

1  1  J  1  o 

<$<N> 

|sxsxsxs?|$(N)> 

o 

1  1  j  k  1  1  O 

(2.19) 


(2.20) 
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Figure  2.5  -  Ground  state  energy/spin  (-E  /NJ)  for  the 

spin  h  XY  model  on  the  squar§  lattice  as  a 
function  of  1/N. 
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and 


<$> 


(N) 


SXSXS^S^ |$(N) > 
1  j  k  1  1  o 


(2.21) 


where  the  sites  i,j,k  and  1  may  be  chosen  from  any  on  the 
N-site  cell.  By  translational  invariance,  the  correlations 
depend  only  on  the  relative  distance  between  these  sites  and 
not  on  their  absolute  location  within  the  cell.  For  two  spin 
correlations,  this  implies  that  there  are  only  a  few  which 
are  independent:  nearest  neighbour,  second  nearest  neighbour, . 
etc.  The  number  of  independent  four  spin  correlations  is 
greater  and  not  as  easily  classified. 

To  evaluate  correlations  numerically,  the  operators 

xy  + 

Si'  i  are  exPressec^  in  terms  of  ladder  operators  through 

the  definitions  (2.6).  The  expressions  are  simplified  by 

(N) 

noting  that  the  ground  state  has  mz=0,  so  that  the 

expectation  value  of  operators  of  the  form  S+S+,  S.S.S  , 

1  3  1  3  k 

S.S.S+S.,  ...,  etc.  in  this  state  is  zero.  Using  this 

1  ]  K  1 

result,  the  correlations  ( 2 . 19 ) - ( 2 . 21 )  become 

(i,j)|<^N)> 

| (il/ jk)  +  ( ik , j 1)  +  (ij ,kl) | >  , 

|  (il,jk)  +  (ik,  jl)  -  (i j ,kl) | $ ^N) >  , 


j<*(N)l 

4  o 

T7T<*(N) 
16  o 

and 

T<<*(N> 
16  o 

where 

(i,j)  s 

and 

(ij ,kl) 

;ts . 

1  3 


;“s+ 

i  3 


(2.22) 

(2.23) 

(2.24) 

(2.25) 


(2.26) 


' 
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The  strategy  for  computing  the  correlations  now  follows  from 
Equation  (2.17),  with  the  operator  A  identified  with  (i,j) 
in  the  case  of  two  spin  correlations  and  with  (il, jk)  +  (ik, jl) ± 
(ij,kl)  for  the  four  spin  correlations  (2.23)  and  (2.24).  The 
procedure  follows  exactly  that  used  to  calculate  the  matrix 
elements  of  -tf  / J,  with  the  single  exception  that  the  states 

A  i 

which  are  counted  are  not  those  which  differ  from  by  two 

nearest-neighbour  spins,  but  those  which  differ  from  ip  ^  by 

a  pair  (or  quadruplet)  of  spins  in  the  set  £  g . ^Ag . ,  which  is 

i 

computed  at  the  beginning  of  the  program.  A  further  complica¬ 
tion  in  the  evaluation  of  (2.24)  is  that  the  contribution  from 
the  states  (i  j  ,kl)  |  ip  ^  >  must  be  counted  with  a  negative  sign, 
but  this  is  easily  accounted  for.  As  the  matrix  elements 

(N)  |A  I  >  are  computed,  the  correlations  ( 2 . 22 )  -  ( 2 . 24  ) 
p  ot 

are  generated  using  the  expansion: 


<4> 


(N) 


IA 


$<N>> 

o 


=£ 

a 


a2<Y(N) 
a  a 


A 


*(N)>  +  2 
Ya 


a ,  3>a 


a  aQ<'l/0 
a  3  3 


(N) 


a  y 


(N) 


> 


a 


(2.27) 


where  the  a  are  the  coefficients  of  the  ground  state 

a 

eigenvector : 

$  (N)  =  £  a  V  (N)  .  (2.28) 

o  a  a 

a 

The  details  of  the  calculation  of  four  spin  correlations  with 
the  form  (2.21)  (or  (2.24))  can  be  obtained  from  the  program 
GXY4C0R  which  has  been  placed  in  the  Statistical  Physics 
Program  and  Data  Library. 

As  examples  of  the  calculations  described  above,  the 
two  and  four  spin  correlations  which  will  be  required  in 
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Chapter  III  have  been  tabulated  in  Table  2.2.  These  corre¬ 
lations  are  expressed  in  terms  of  Pauli  operators,  cr^, 
which  are  related  to  the  spin  h  operators  by  S^=  (fi=l) , 

so  that 

<axax>  =  4<SXSX>  ,  (2.29) 

1313 

and  <axaxa^a^>  =  16  <SXSXS^S^>  ,  (2.30) 

i  3  k  1  ljkl 

and  where  of  course  the  expectation  value  is  taken  in  the 
ground  state  •  If  6  is  the  lattice  spacing  then  second- 

nearest  neighbour  sites  are  /2S  units  apart  while  third-nearest 
neighbours  are  25  units  distant,  on  the  square  lattice.  In 
Table  2.2(a),  rows  1  and  2  represent  respectively  2nd  and 
3rd  nearest  neighbour  correlations.  The  four  spin  correla¬ 
tions  involve  sites  lying  on  a "domino"  formed  by  adjacent 
unit  squares,  or  plaquettes.  The  exact  geometry  is  not 
important  here  but  will  be  discussed  in  the  next  chapter.  The 
first  two  rows  of  Table  2.2(b)  are  nearest  neighbour  and 
second-nearest  neighbour  correlations  on  the  triangular 
lattice  while  the  four  spin  correlations  are  again  taken 
over  sites  of  a  "domino"  formed  by  two  adjacent  triangular 
plaquettes.  Again,  this  will  be  made  explicit  in  Chapter  III. 
The  two  spin  correlations  were  calculated  by  Dr.  Leo  Marland 
and  obtained  from  him  {2.6},  as  was  the  data,  e.g.  basis 
states,  ground  state  eigenvector,  required  to  calculate  the 
four  spin  correlations.  It  should  be  mentioned  here,  if  it 
is  not  already  apparent, that  all  the  programs  described  above 
will  work  without  modification  for  cells  on  the  triangular 
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Table  2.2  -  Ground  state  correlations  for  the  spin  h  XY  model 

on  the  (a)  square  and  (b)  triangular  lattices. 


(a) 


N  8  10  16  18 


x  X  . 

0°  m 

0.52273 

0.50746 

0.48685 

0.48797 

xa*  > 
o  26 

0.52273 

0.50746 

0.48685 

0.47278 

0.09091 

0.09919 

0.11278 

0.08842 

AAA> 

0.09091 

0.07222 

0.04843 

0.05299 

x  x  v  y^ 
lQ3a2a6 

0.09091 

0.09919 

0.07997 

0.08264 

^  x  y  y^ 
Ia3a4a6 

0.09091 

0.09919 

0.07997 

0.06909 

(b) 


13 


19 


N 


1.  <0 

2 .  <o 

3 .  <o 

4 .  <a 


x  x. 

ax> 
o  o 


xx  . 
o  /36 


xx  xx. 
la2a3°4> 


x  x  y  V 

2a4aIa!> 


0.57143 


0.56103 

0.52111 

0.47837 


0.55052  0.54391 
0.50295  0.48914 
0.45715 


0.13100 


0.11796 
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lattice.  All  that  is  needed  to  prime  the  sequence  is  a 
choice  of  finite  cell  and  the  elements  of  the  symmetry  groups 
Tn  and  PN  expressed  as  permutations . 

Figures  2.6  and  2.7  are  respectively  graphs  of  the 
two  spin  and  four  spin  correlations,  for  both  lattices,  as 
functions  of  1/N.  The  two  spin  correlations  on  the  triangu¬ 
lar  lattice  are  well  behaved,  with  the  second-nearest  neighbour 
correlation  forming  a  near  perfect  straight  line.  There  is 
a  slight  deviation  of  the  nearest  neighbour  correlation  at 
N=19  which,  since  the  ground  state  energy  is  proportional  to 
the  nearest-neighbour  correlation,  is  also  evident  in  the 
former  quantity  and  is  of  unknown  origin  {1.17}.  A  rather 
surprising  feature  of  the  square  lattice  data  is  that  both 
second  and  third-nearest  neighbour  correlations  are  equal 
for  N£l6.  This  degeneracy  is  absent  in  the  results  from  the 
18-spin  cell  calculation,  which  is  to  be  expected  as  these 
correlations  are  certainly  not  the  same  in  the  infinite 
lattice  limit.  This  phenomenon  may  be  due  to  some  unaccounted 
symmetry  of  the  N£l6  cells,  but  whatever  its  cause,  it  is 
clear  that  extrapolation  to  the  infinite  lattice  will  be 
inaccurate  if  this  is  based  on  either  all  data  points  or 
simply  the  results  for  N^16.  When  the  four  spin  correlations 
are  examined,  the  situation  appears  much  worse,  for  the  square 
lattice  at  least.  Degeneracy  is  evident  again:  for  N=8,  all 
four  of  the  four  spin  correlations  are  equal,  for  N=10  three 
of  them  are  the  same  while  at  N=16  two  of  them  are  the  same, 
and  it  is  not  until  N=18  that  the  four  correlations  are  all 
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Figure  2.6  -  Selected  two  spin  correlations  for  the  spin  h 

XY  model  on  the  square  (squares)  and  triangular 
(triangles)  lattices  as  functions  of  1/N. 
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Figure  2.7  -  Selected  four  spin  correlations  for  the  spin  ^ 

XY  model  on  the  square  (squares)  and  triangular 
(triangles)  lattices  as  functions  of  1/N;  (see 
Figures  3.5  and  3.6  for  a  definition  of  the 
coordinates ) . 
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distinct.  Aside  from  this  however,  and  perhaps  acting  as  a 
cause  for  some  of  the  degeneracy,  is  the  fact  that  for  N£l0, 
the  area  of  the  cell  is  of  the  same  order  as  the  area  over 
which  the  four  site  correlations  are  evaluated.  For  the  same 
reason  that  the  ground  state  energy  per  spin  calculated  for 
N=4  (and  5)  does  not  lie  on  the  straight  line  formed  by  the 
larger  N  results,  so  to  it  is  expected  that  the  four  spin 
correlations  for  N=8  and  10  will  likely  not  lie  on  a  line  with 
larger  N  data.  All  the  spins  of  cells  with  a  small  number  of 
sites  are  too  highly  correlated,  because  of  the  periodic 
boundary  conditions,  to  accurately  reflect  the  situation  in 
the  infinite  lattice  over  distances  greater  than  ^  half  the 
cell's  linear  dimension.  With  this  criterion  the  four  spin 
correlation  for  the  16  spin  cell  is  just  within  the  limit  of 
applicability,  although  too  much  trust  should  not  be  placed 
in  results  for  N>16.  For  example  both  the  correlations 

and  <a^a^a-|a^>  extrapolate  to  negative  numbers 
using  only  the  data  points  N=16  and  N=18!  Although  only  two 
data  points  are  available  for  the  four  spin  correlations  on 
the  triangular  lattice,  there  is  reason  to  believe  that  an 
extropolation  to  the  infinite  lattice  based  on  these  two 
points  alone  would  be  sufficiently  accurate.  First,  the 
two  spin  correlations,  particularly  the  second-nearest 
neighbour  correlation,  are  well  behaved,  showing  no  degeneracy 
which  would  "artificially"  shift  some  data  points.  Second, 
the  area  over  which  the  four  site  correlation  is  taken  is  not 
large  with  respect  to  cell  size,  even  at  N=9,  for  which  the 
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ratio  of  these  two  areas  is  roughly  comparable  with  that  for 
N=16  on  the  square  lattice.  The  maximum  distance  between 
points  on  one  triangular  "domino"  is  /36,  compared  with 
a  maximum  separation  of  /5/2S  for  points  of  the  "domino"  on 
the  square  lattice.  Although  the  four  spin  correlations  for 
N=19  should  be  calculated,  it  is  to  be  expected  that  they  will 
lie  close  to  the  lines  formed  by  the  data  points  already 
known  and  plotted  in  Figure  2.7. 


- 


CHAPTER  III 


VORTICES  IN  THE  GROUND  STATE  OF  THE  SPIN  h  XY  MODEL 
A.  Vortices 

Let  u(x,y,z)  be  a  vector  field,  or  flow,  defined  on 
a  three-dimensional  Cartesian  space.  The  vorticity  of  this 
flow,  w,  is  defined  by  w=curl  u.  The  usefulness  of  this 
quantity  in  understanding  flow  configurations  derives  in 
great  part  from  the  fact  that  vorticity  corresponds  to  a 
rotation  of  the  vector  field.  This  follows  from  the 
definition  of  curl  as 

h*curl  u  =  lim  u*dl  ,  (3.1) 

s-*0 

where  n  is  the  unit  vector  normal  to  the  surface  S,  around 
the  edge  of  which  the  line  integral  is  taken.  In  contrast  with 
this  it  is  possible  to  make  the  definition  adopted  in  this  work: 
a  configuration  is  said  to  possess  unit  vorticity  (antivortici- 
ty)  about  a  selected  point  if  on  traversing  a  contour  about 
that  point  in  a  positive  direction  the  vector  field  makes  a 
complete  ( 2 tt  rad.)  turn  in  the  same  (opposite)  direction.  It 
is  clear  that  this  definition  is  appropriate  only  for  fields 
restricted  to  a  plane.  Illustrations  of  this  definition  are 
presented  in  Figure  3.1.  The  set  of  flow  configurations  which 
possess  vorticity  m  is  not  the  same  as  the  set  of  configurations 
which  possess  vorticity  according  to  the  second  definition.  It 
is  possible  to  exhibit  flows  which  possess  vorticity  according 
to  either  definition,  both  definitions  or,  of  course,  which  do 
not  have  any  vorticity  at  all.  For  example,  it  is  clear  that 
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the  first  of  Figure  3.1a)  is  such  that  cu=  0  at  every  point 
while  in  general  for  the  second  configuration.  This 

illustrates  one  of  the  most  useful  properties  of  the  opera¬ 
tional  definition ,  namely  that  a  vortex  (or  antivortex)  remains 
so  if  the  vector  field  u(x,y)  is  subject  to  rotation  by  a  fixed 
amount  at  each  point  (x,y)  .  The  vorticity  oj  does  not  share 
this  invariance. 

Consider  now  the  plane  rotator  model,  defined  by  the 
Hamiltonian 

H  =  -J  Z  S.-S.  =  -j  z  cos  (tj) .-(}). )  (3.2) 

rK  ....in  ....  in 

<i,3>  J  <i,3>  J 

where  H  S  Jl  =1  and  <Jk  is  the  angle  makes  with  some  arbitrary 
axis,  0£4k£2tt  .  This  model  is  the  classical  analogue  of  the 
spin  h  XY  model  and  is  in  fact  the  same  obtained  from  the 
spin  s  XY  model  in  the  limit  as  s->°°.  To  help  illustrate  the 
concept  of  vorticity  before  proceeding  to  a  discussion  of 
vortices  in  the  spin  h  XY  model,  it  will  be  shown  below  that 
the  probability  of  finding  a  vortex  (or  an  antivortex)  about 
any  given  location  in  the  square  lattice,  in  the  limit  of 
infinite  temperature,  is  0.167*  It  should  be  mentioned  here 
that  consideration  of  a  vector  field,  S^,  defined  on  a  lattice 
does  not  interfere  with  the  definition  of  vorticity  adopted 
above.  The  effect  of  the  lattice  is  to  introduce  a  contour 


44 


\\//  ^  ^  \ 

^  \  /  /  \  \ 

\  \  /  / 

/  /  \  \  \  ^  / 


(a)  unit  vortices 


\  ^  ^  / 

\  \  /  / 
1  /  \  \ 


(b)  unit  antivortices 


Figure  3.1  -  Vortex  configurations  in  a  two-dimensional  vector 

field. 
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of  minimum  size  which,  on  the  square  lattice,  is  just  a 
square  with  an  edge  of  length  6  and  vertices  formed  by  lattice 
sites.  It  is  about  this  contour,  known  as  a  plaquette,  that 
vorticity  will  be  calculated. 

In  the  limit  of  infinite  temperature,  the  probability 
of  finding  the  angle  of  the  spin  S^,  with  respect  to  some 
fixed  axis,  lying  in  the  range  0  to  9  +  d0  is  just  d@/27T, 
independent  of  the  orientation  of  any  other  spin  on  the  lattice. 
The  problem  of  finding  the  probability  of  a  vortex  configura¬ 
tion  then  reduces  to  that  of  finding  the  proportion  of  all 
spin  configurations  on  a  plaquette  which  are  vortices. 

Further,  because  a  vortex  is  invariant  under  a  rotation  of  all 
spins  by  a  fixed  amount,  it  is  only  necessary  to  consider  the 
orientation  of  three  spins  relative  to  that  of  the  fourth, 
which  can  be  assumed  fixed.  This  fixed  spin  will  be  placed 
on  site  1,  and  the  positive  direction  along  the  plaquette 
will  be  chosen  to  coincide  with  the  direction  1234.  The 
spin  on  site  3  will  be  placed  at  angle  of  9  with  respect  to 
the  axis  defined  by  the  spin  on  site  1.  Two  cases  will  be 
distinguished:  O£0£tt,  illustrated  by  Figure  3.2(a)  and 
tt£0<2tt,  illustrated  by  Figure  3.2(b).  Once  the  spins  at 
sites  1  and  3  are  fixed,  the  question  becomes  the  determination 
of  the  orientations  of  the  spins  on  sites  2  and  4  which  will 
yield  a  vortex  configuration.  Consider  first  Figure  3.2(a). 

The  spin  on  site  2  may  lie  anywhere  in  the  hatched  region, 
including  the  axis,  while  the  spin  on  site  4  may  lie  in  the 
hatched  region  about  this  site,  that  is  at  an  angle  tt < < tth- 0 
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Figure  3.2  -  Calculation  of  the  vorticity  of  the  plane  rotator 

model,  in  the  infinite  temperature  limit,  on  the 
square  (a),  (b)  and  triangular  (c)  lattices. 
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with  respect  to  the  axis  defined  by  spin  1.  These  hatched 
regions  are  constructed  so  that  no  spin  lying  inside  these 
regions  is  more  than  +tt  rad.  "away"  from  the  axes  defined  by 
both  the  spins  on  site  1  and  site  3.  This  ensures  the 
spins  make  a  positive  turn  on  going  around  the  plaquette 
in  a  positive  direction.  Note  that  the  results  imply  that  no 
vortex  can  be  formed  when  both  spin  1  and  spin  3  are  parallel 
(0=0) .  The  proportion  of  all  spin  configurations  which  are 
vortices  when  spin  3  lies  in  the  range  0  to  0+d0,  with 
0^0^71,  is  then 


dP 


(a) 


0 _ 

2tt 


sites:  4 


(2tt-0) 

2tt 

i _ i 

2 


(3.3) 


Similarity,  the  proportion  of  all  configurations  which  are 
vortices  when  spin  3  lies  in  the  range  0  to  0+d0,  with 
7t£0s?2it,  is,  from  Figure  3.2(b), 


dP 


(b) 


( 2 7T—  0  )  d0  Q_ 


2  7T 
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2tt  2tt 

3  4 


(3.4) 


It  follows  from  previous  remarks  that  the  probability  of 
finding  a  vortex  on  a  given  plaquette,  in  the  infinite 
temperature  limit  is 
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P 
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2  TT 

+  /  dP 


(b) 


TT 

/  e(27T-0)d0 
o 


=  0.167  . 


(3.5) 


The  probability  of  finding  an  antivortex  on  a  given  plaquette 
is  also  found  to  be  0.167,  a  result  which  follows  immediately 
from  the  above  arguments  once  it  is  noted  that  an  antivortex 
can  always  be  formed  from  a  vortex  by  inverting  two  diagonally 
opposite  spins,  e.g.,  the  spins  on  sites  1  and  3.  The 
probability  of  finding  either  a  vortex  or  an  antivortex  in 
the  infinite  temperature  limit  of  the  plane  rotator  model  is 
then  0 . 333. 

The  same  type  of  calculation  can  be  carried  out  for 
the  plane  rotator  on  the  triangular  lattice,  for  which  the 
basic  plaquette  is  a  triangle  of  side  6  with  vertices  formed 
by  lattice  sites.  Figure  3.2(c)  illustrates  the  situation 
in  this  case.  Positive  direction  around  the  plaquette  is 
taken  as  1  2  3 .  The  fixed  spin  is  placed  on  site  1  and  the 
spin  on  site  2  is  oriented  at  an  angle  of  0,  O£0£tt,  with 
respect  to  the  axis  defined  by  the  fixed  spin.  Note  that 
for  the  triangular  lattice  vortex  configurations  correspond 
to  O^0£tt  while  antivortex  configurations  are  formed  for 
7T£0<2tt.  Referring  to  Figure  3.2(c),  the  differential 
probability  of  finding  a  vortex  on  a  given  plaquette  is 
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so  that  the  total  probability  is 

77  (  \ 

P  =  /  dP  C 
v  J 
o 

1  77 

=  j  0d6 
4tt  o 

=  0.125  ,  (3.7) 

as  in  the  case  of  the  square  lattice.  The  probability  of 
finding  an  antivortex  is  then  0.125,  so  that  the  total 
probability  of  finding  a  vortex  or  an  antivortex  on  the 
triangular  lattice  is  0.250. 

B .  Construction  of  the  Vortex  Operators 

In  the  previous  section  the  concept  of  vorticity  was 
introduced  and,  almost  immediately,  adopted  from  vector 
fields  defined  on  a  continuous  space  to  those  defined  on  a 
two-dimensional  lattice.  One  further  modification  must  be 
introduced  before  vortices  in  the  ground  state  of  the  spin  h 
XY  model  can  be  discussed,  and  this  is  the  quantization  of 
the  vector  field,  S..  The  considerations  outlined  above 

l 

were  appropriate  for  a  classical  field  which  could  assume 
any  in  a  continuous  range  of  orientations.  In  the  spin  h  XY 
model,  S.  becomes  a  quantum  mechanical,  spin  h  operator;  what 
becomes  of  the  vorticity? 

The  answer  to  this  question  adopted  in  this  work  is  to 
construct  an  operator,  the  vorticity  operator  v,  the 
eigenstates  of  which  are  a  vortex  or  an  antivortex  configura¬ 
tion,  or  neither,  over  a  basic  plaquette.  The  corresponding 
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eigenvalues  will  be  +1,  -1  and  0  respectively.  The 
expectation  value  of  this  operator  in  the  ground  state  of 
the  spin  h  XY  model  will  then  be  interpreted  as  representing 
the  mean  number  of  vortices  per  site  present  in  this  state. 

To  construct  the  vorticity  operator  appropriate  for 
the  square  lattice,  consider  first  the  basic  plaquette  shown 
in  Figure  3.3.  The  spins  on  sites  1  and  3  are  quantized 
along  the  x  axis  while  the  spins  on  sites  2  and  4  are  quantized 
along  the  y  axis.  The  individual  site  operators  are  Pauli 


operators  a,  and  so  there  are  a  total  of  16  distinct  states 

possible  on  each  plaquette.  Since  the  vorticity  is  invariant 

under  a  spin  flip  applied  to  all  sites,  only  8  of  these 

states  will  contribute  independently  to  the  construction  of  v. 

These  8  states  are  shown  in  Figure  3.3.  Of  these  8  states, 

only  state  1  represents  a  vortex  while  only  state  7  represents 

an  antivortex;  all  others  show  no  vorticity  and  so  v  applied 

to  these  states  must  be  zero.  The  next  step  is  to  pick  out 

X  y  x  v 

eight  functions  of  the  four  spin  variables  0  2'0  2'  and 

from  which  to  construct  an  operator  v  which  gives  +1  when 
applied  to  state  1,  -1  when  applied  to  state  7  and  annihilates 
all  the  other  states.  These  functions  must  be  even  functions 
of  the  basic  operators,  a  suitable  choice  being  o^a^o^a^, 

a2°3'  a3°4'  ala3'  a2a4'  al°2a3a4  an<^  T*ie  fanal  task  is 
to  solve  the  8  equations  for  the  coefficients  of  these 

operators  in  the  expansion  of  v  obtained  from  each  of  the  8 

states  in  Figure  3.3.  Passing  over  the  details  of  this  step. 


it  is  found  that  only  four  of  the  above  operators  are  needed 


51 


/. 


Cv7  -  I 


2. 


V 


h. 


<^o>  '0 


<u>  •  0 


«0 


S\ 


6. 


7. 


g. 


<u> 


--  O 


if 


•  o 


=  ") 


<u>  =0 


Figure  3.3  -  Vortex  states  on  the  square  lattice. 
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to  describe  v,  which  can  then  be  written  as 


1 ,  x  y  x  y 


v  =  4(°i°2  '  °la4  +  a304 


x  y  y  x 

rr-4  —  rr-4  n 


a2a3> 


(3.8) 


When  the  expectation  of  this  operator  is  taken  for  the  spin  h 
XY  model  however,  it  is  found  that  <v>=0,  at  any  temperature, 
since  <a^cn>=0  by  symmetry,  for  any  sites  i,j  on  the  lattice. 
The  interpretation  of  this  result  is  that  in  the  ground  state 
the  number  of  vortices  and  anti vortices  is  equal.  To  circum¬ 
vent  this  result  the  operator  v  is  squared  and  the  subsequent 
expectation  value  interpreted  as  the  (mean)  total  number  of 
vortices  and  anti vortices  per  site  in  the  ground  state. 
Carrying  out  this  operation  yields 


2  1  , yy  xx,  x  y  x  y . 

V  =  4(1  ~  °2°4  "  al°3  al°2a3°4) 


(3.9) 


so  that 

^  2^  1,,  ~xx  ^  ^  x  x  y  y,.  , 

<v  >  =  j(l  -  2<a0a/25>  +  ala3a2a4  ^  '  (3.10) 

where  the  equality  <aia3>=<a2a4>  holds,  at  any  temperature,  in 
the  XY  model  due  to  the  symmetry  of  x  and  y  in  and 

<0^0^02 a4>=<aia2a3a4>  s;*-nce  t^ie  operators  at  different  sites 
commute . 

The  operator  v1  appropriate  for  the  triangular  lattice 
can  be  constructed  in  a  similar  manner.  In  Figure  3.4(a)  the 
possible  states  on  a  triangular  plaquette  are  illustrated. 

The  spin-flip  states  are  again  not  necessary.  On  this  lattice 
it  is  necessary  to  choose  the  direction  of  quantization  of  o ^ 
to  be  the  v-direction  and  that  of  to  be  the  u-direction. 


s 
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Figure  3.4  -  Vortex  states 


on 


the  triangular  lattice. 
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Neither  of  these  are  orthogonal  to  x  or  y,  and  in  fact  the 
quantities  GU  and  gV  can  be  written  in  components  as 

ou  =  /T/2  oy  -  j  ox  ,  (3.11a) 

ov  =  -/3/2  0y  -  j  ox  .  (3.11b) 

The  first  state,  for  which  <v'>=l,  has  all  three  spins  oriented 
in  the  positive  coordinate  directions  while  the  other  three 
have  at  least  one  spin  oriented  in  the  negative  coordinate 
direction.  The  operator  v*  can  thus  be  written  down  immediately 
as 


,  1,.  ,  x  u  .  x  v  ,  v  u.  /0  . 

v'  =  j(l  +  g1g3  +  0^2  +  °2°3^  *  (3.12) 

2 

This  operator  has  the  interesting  properties  that  v,=(v')  and 
< v '  > 7^0 .  The  former  can  be  shown  by  direct  computation  while 
the  latter  can  be  seen  in  the  following  manner.  Using  the 
relations  (3.9),  the  identities  <GXo^>=0,  <oXoX>=<oY oY >  and 
translational  invariance,  the  ground  state  expectation  of 
(3.12)  is 


< v 1  >  =  +  <ai  (”^3)  >+<ai  (”^3)  >+<  (J_fa2)  ( 


2^3 


)> 


1  3  x  x 

T  “  7T<a  G  r  > 

4  806 


(3.13) 


To  understand  why  this  quantity  is  non-zero  one  must  consider 
the  number  of  vortices  and  antivortices  that  are  present  in  a 
lattice  tiled  with  plaquettes  of  the  form  in  Figure  3.4(a). 
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But  first  it  is  necessary  to  generate  antivortices,  which  have 
not  yet  been  obtained  on  the  triangular  lattice,  the  operator 
v'  having  only  +1  and  0  as  eigenvalues.  The  accomplishment 
of  this  task  is  not  difficult.  It  is  in  fact  not  possible  to 
construct  a  triangular  lattice  in  which  all  the  plaquettes 
will  have  axes  of  quantization  as  shown  in  Figure  3.4(a).  The 
plaquettes  which  form  "in-between"  these  on  the  bulk  lattice 
will  have  axes  of  quantization  as  shown  in  Figure  3.4(b), 
after  a  rotation  of  each  spin  to  bring  that  on  site  1  into 
the  orientation  shown.  For  clarification  of  the  geometry 
involved,  see  Appendix  III,  where  it  is  shown  that  the  ground 
state  of  the  antiferromagnetic  plane  rotator  on  the  triangular 
lattice  consists  of  closely  packed  vortices  and  antivortices 
of  the  form  shown  in  Figure  3.4(a)  1  and  (b)  1.  From 
arguments  identical  to  that  which  lead  to  (3.12),  the  anti¬ 
vortex  operator  v"  is  found  to  be 

v"  =  —4(1  +  a]_a3  +  °i°2  +  a2a3  *  (3.14) 

2 

Notice  that  again  v"=(v")  and  <v">^0.  In  fact  <v">=<v'>,  so 

that  the  interpretation  of  this  operator  expectation  value  is 

that  it  gives  the  (mean)  total  number  of  vortices  and  anti- 

2 

vortices  in  the  ground  state  per  plaquette,  just  as  <v  >  does 
on  the  square  lattice.  In  addition  there  is  an  extra  factor 
which  arises  from  the  geometry  of  the  triangular  lattice. 
Unlike  the  situation  on  the  square  lattice  where  there  is  one 
square  plaquette  per  lattice  site,  there  are  two  plaquettes 
per  lattice  site  on  the  triangular  lattice.  The  total  number 
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of  vortices  and  antivortices  in  the  ground  state  per  lattice 
site  is  then  given  by 

<V’>  =  I  -  htt  -  (3-l5) 

and  it  is  this  expectation  value  which  will  be  exhibited  as 
a  function  of  the  number  of  spins  per  cell,  N,  in  the  sequel. 
This  expression  involves  only  the  nearest-neighbour  spin 
correlation  and  so  is  proportional  to  the  ground  state  energy 
per  spin. 

C .  Construction  of  the  Vortex-Antivortex  Pair  Operators 

Now  that  the  vortex  operators  applicable  to  the  square 
and  triangular  lattices  have  been  derived,  attention  can  be 
focused  on  other  operators  that  may  be  interesting,  in 
particular  the  vortex-antivortex  pair  operator.  The  eigen¬ 
states  of  this  operator,  with  non-zero  eigenvalue,  consist  of 
a  vortex  configuration  about  one  plaquette  and  an  antivortex 
about  the  adjacent  plaquette.  All  other  states,  including 
those  consisting  of  adjacent  vortices  or  antivortices,  are 
assigned  zero  eigenvalue.  This  operator,  or  rather  the 
expectation  value  of  this  operator,  is  useful  for  the 
estimation  of  the  number  of  free,  unbound  vortices  in  the 
ground  state  of  the  spin  h  XY  model. 

The  pair  operator  on  the  triangular  lattice  will  be 
considered  first,  since  it  is  simpler  to  derive  than  the 
corresponding  operator  on  the  square  lattice.  The  eight 
configurations  which  must  be  considered  to  derive  this 
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operator  are  shown  in  Figure  3.5.  The  contour  over  which 
these  configurations  are  defined  is  formed  by  joining  two 
adjacent  plaquettes  and  will  be  termed  a  domino  both  for  the 
triangular  and  square  lattices.  Unlike  the  situation  which 
obtained  when  the  vortex  operator  was  constructed  on  the 
triangular  lattice,  there  is  no  need  to  consider  any  other 
possible  domino  configurations  since  all  such  configurations 
can  be  made  similar  to  those  shown  in  Figure  3.5  by  transfor¬ 
mations  in  the  space  group  of  the  lattice  and/or  a  uniform 
rotation  of  all  the  spins.  The  pair  operator  p'  can  then  be 
represented  as 


i  _ 


=  t  (1  + 


XU,  u  V  ,  V  X  , 

°la2  +  °2°3  +  a3al  + 


XU,  u  u 
al°4  +  <J2°4  + 


V  U  ,  X  U  V  u, 

o,o,  +  ala2a3a4 


3  4 


(3.16) 


This  formula  can  be  justified  by  the  following  argument  in 
which  all  numbers  refer  to  Figure  3.5,  a=x,u,v  and  i=l,2,3,4. 

In  the  configuration  1,  all  g^  are  positive  so  <p ' >  is  1. 

When  one  g^  is  negative,  as  in  configurations  2,  3  and  5,  the 
four  spin  term  and  three  of  the  two  spin  terms  in  p'  will  be 
-1,  while  the  three  remaining  two  spin  terms  are  +1  which  gives 
<p'>=0.  When  two  of  the  g^  are  negative,  as  in  configurations 
4,  6  and  7,  four  of  the  six  spin  terms  (those  containing  at 
least  one  of  the  negative  a^)  will  be  -1,  while  the  other  terms 
in  p'  will  be  +1  so  that  <p'>  is  again  zero.  Finally  when  three 
of  the  g^  are  negative,  as  in  configuration  8,  the  four  spin 
term  is  negative  as  is  each  of  the  three  two  spin  terms  con¬ 
taining  the  one  remaining  positive  ;  all  other  terms  will 
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Figure  3.5  -  Vortex-antivortex  pair  states  on  the  triangular 

lattice.  (vEvortex,  aEantivortex) 
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again  be  +1  so  that  <p'>=0  in  this  case  as  well.  It  is  also 

2 

straightforward  to  observe  that  p'=(p')  .  The  implication  of 

this  result  for  the  interpretation  of  the  expectation  value 
of  p  is  that  p  does  not  discriminate  between  a  pair  which  is  a 
v(ortex)-a  (ntivortex)  and  one  which  is  a.-v.  This  result 
is  in  fact  obvious  from  the  derivation  of  p*.  In  calculating 
<p*>,  use  is  made  of  the  following  relations: 


.  X  U  ,  U  V  ,  V  x^ 

<0l°2+CT2a3+a3al>  = 


3^  x  x^ 
-Tr<a  af> 
2  o  6 


(3.17) 


XU  V  u 
<ala4>=<a3a4>  = 


1  x  x^ 
-7r<a  ax> 

2  o  o 


(3.18) 


u  u 
<a2a4> 


X  X  ^ 

=  <0o?36> 


(3.19) 


In  addition,  the  expectation  value  <aga2a3a4>  can  simplified 
before  evaluation  by  rotating  the  axes  u,  v,  x  by  120° 
according  to  x-*v,  v-*u  and  u-*x.  The  ground  state  expectation 
value  is  invariant  under  this  transformation  since  is 

invariant  under  this,  and  any,  rotation  in  the  x-y,  or  x-u-v, 
plane;  the  desired  four  spin  correlation  becomes 


^  v  x  u  x^  1^  x  x  x  x^  3^  x  x  y  y^ 
aia2a3a4  ~  4  aia2a3a4  4  a2a4ai°3 


(3.20) 


Using  the  results  ( 3 . 17 ) -  (  3 . 20 )  ,  the  expectation  value  of  p* 


can  be  written  as 


.  1  5  xx  ,  1  x  x  ^.l^xxxx^  3^.xxyy^ 

<p  >=8  -  16  aoCT6  +  8<°oa/36>  +  32<0ia2a3a4>  ~  32<a2a4°la3> 


(3.21) 
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This  result  measures  the  (mean)  number  of  vortex-antivortex 
pairs  per  domino  ,  in  the  ground  state  of  the  spin  %  XY  model 
on  the  triangular  lattice.  On  this  lattice  however,  there  are 
three  different  domino  patterns  per  lattice  site  so  that  the 
expectation  value  which  will  be  reported  in  the  next  section, 
giving  the  number  of  vortex-antivortex  pairs  per  lattice  site  / 
is 


(3.22) 


The  pair  operator  p  which  is  applicable  to  the  square 


lattice  can  of  course  be  derived  from  similar  considerations. 

In  this  case  however  the  domino  made  up  of  two  adjacent 

square  plaquettes  contains  six  lattice  sites  and  even  though 

the  spin  on  one  site  may  be  fixed  there  are  still  32  pair 

states  which  must  be  considered  in  the  derivation  of  p. 

Several  of  these  are  shown  in  Figure  3.6.  It  is  actually  an 

2 

easier  task  to  derive  p  directly  rather  than  construct  p  and 
square  the  resulting  operator.  At  any  rate  the  procedure  is 
straightforward,  if  somewhat  tedious  and  the  final  result  can 
be  written  down  without  consideration  of  the  details  of 
derivation  as 


x  x  y  v  _  x  x  v  y 


a3a4Q6  ° la5a2a  4 


a, a-ata 


(3.23) 
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Figure  3.6  -  Selected  vortex-antivortex  pair  states  on  the 

square  lattice.  (vEvortex,  aEantivortex) 
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where  the  numbering  is  given  in  Figure  3.6.  In  taking  the 
expectation  value  of  (3.23)  notice  must  be  taken  of  the  number 
of  dominoes  per  lattice  site,  on  an  N  site  cell,  which  on  the 
square  lattice  is  just  two.  Including  this  correction  factor 
leads  to  the  expectation  value 


x  x  y  y 


V  V  V  V 


(3.24) 


+  2<a,a0aia^> 
13  2  6 


which  represents  the  (mean)  number  of  vortex-  antivortex 
pairs,  per  lattice  site,  in  the  ground  state  of  the  spin  h  XY 
model  on  the  square  lattice. 

D .  Presentation  and  Discussion  of  Results 

It  is  a  simple  matter  now  to  match  the  expressions  for 
2  2 

<v  >,  <v'>,  <p'>  and  <p  >  given  above  and  the  various  two  and 
four  spin  correlations  presented  in  Chapter  2,  section  E,  to 
determine  the  vorticity  and  pair  vorticity  operators  as  func¬ 
tions  of  1/N.  These  are  presented  in  Table  3.1  and,  graphically, 
in  Figures  3.7  and  3.8.  Notice  that  considerations  of  time 
and  expense  have  limited  the  calculation  of  <p'>  to  the  9  and 
13  spin  cells.  This  quantity  is  not  presented  for  the  7  spin 
cell  because  this  particular  cell  does  not  allow  the  calcula¬ 
tion  of  second-nearest  neighbour  correlations,  every  site  being 

a  nearest  neighbour  to  each  other  site. 

2 

Both  <v'>  and  <v  >  are  seen  to  behave  in  an  approximately 
linear  fashion  as  functions  of  1/N,  although  there  is  a 
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Table  3.1  -  Vorticity  and  pair  vorticity  for  the  spin  h  XY  model 

on  the  (a)  square  and  (b)  triangular  lattices. 


(a) 


N 


8 


10 


16 


18 


<v2>  0.01134 

<p2>xl03  5.68 


v 


free 


^0 


0.01432  0.01861 

3.79  9.46 

3.37xl0“3  ^0 


0.01926 

9.62 

%0 


(b) 


N 


7 


13 


19 


A 

> 

V 

0.07144 

0.07922 

0.08710 

0.09206 

A 

V 

— 

0.05244 

0.05718 

— 

vf  ree 

- 

-0.01283 

-0.01363 

- 

— 

0.00691 

0.00774 

— 

<v2> 
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0.122 
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0.082 


0.062 


0.042 


0.022 


1/N 

Figure  3.7  -  Expectation  values  of  vortex  operators  for  the 

square  (squares,  left-hand  scale)  and  triangular 
(triangles,  right-hand  scale)  as  functions  of  1/N. 
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Figure  3.8  -  Expectation  values  of  vortex-antivortex  pair 

operators  for  the  square  (squares,  left-hand 
scale)  and  triangular  (triangles,  right-hand 
scale)  lattices  as  functions  of  1/N. 
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distinct  curvature  to  the  results,  especially  noticeable  for 

<v*>.  Because  of  this  the  best  fit  line  used  to  extrapolate 

these  quantities  was  drawn  with  a  bias  in  favor  of  the  larger 

N  data  and  the  estimated  uncertainty  in  these  extrapolations 

2 

raised  accordingly.  For  <v  >  the  value  found  for  the  infinite 

square  lattice  is  0.025±0.002,  while  the  value  of  <v' >  for  the 

infinite  triangular  lattice  is  found  to  be  0.104+0.008.  It 

seems  at  this  point  to  be  coincidental  that  the  latter  number 

is  slightly  more  than  four  times  the  former  extrapolation. 

2 

The  situation  with  regard  to  <p  >  and  <p'>  is  less  certain. 

Only  two  data  points  are  available  from  which  to  determine  a 
value  of  <p'>  in  the  infinite  lattice  limit.  The  value 
obtained  for  this  quantity,  0 . 067±0 . 010,  has  a  correspondingly 
larger  uncertainty.  The  data  varies  widely  but  an  extrapola¬ 
tion  of  <p2>  to  the  infinite  lattice  case  of  (1. 23±0 . 10) xlO  ^ 

can  be  quoted,  with  very  limited  accuracy.  With  regard  to  this 

2 

curve  however  it  is  interesting  to  note  that  if  <p  were 

doubled,  the  improvement  in  the  data  would  be  such  as  to  allow 

an  extrapolation  to  be  made  with  some  confidence,  although 

the  same  negative  curvature  of  the  data  as  observed  in  the 

vorticity  operators  would  be  present.  Needless  to  say  the 

2 

calculations  leading  to  <p  >^q  have  been  checked  many  times  for 
accuracy.  If  some  error  is  present,  it  is  likely  due  to  a 
geometric  factor,  a  result  of  the  fact  that  of  the  four  even  N 
cells  on  the  square  lattice,  only  the  N=10  cell  was  not 
reflection  invariant  and  so  has  a  smaller  symmetry  group  than 
the  others.  More  will  be  said  on  this  subsequently. 
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The  interpretation  of  the  results  presented  above 

begins  with  examination  of  the  infinite  temperature  limit  of 

2  2 

the  operator  expectation  values  <v  >,  <v*>,  <p  >  and  <p’>  . 

In  this  limit  all  correlations  are  zero  and  the  components  of 
the  site  spin  operators  are  distributed  randomly.  The 
vorticity  per  lattice  site  on  the  square  lattice  is  then  0.250. 
This  in  turn  means  that  the  number  of  vortices  per  lattice 
site,  or  equivalently,  per  plaquette  is  0.125,  which  is  not  in 
agreement  with  the  infinite  temperature  limit  obtained  in 
section  A  for  the  classical  (s=°°)  XY  model.  The  expectation 
value  of  v* ,  in  this  limit,  is  0.500,  which  has  an  interpreta¬ 
tion  as  the  average  number  of  vortices  and  antivortices,  per 
lattice  site,  on  the  triangular  lattice.  Recalling  that  on 
this  lattice  there  are  2  plaquettes  per  site  and  that  the 
number  of  vortices  is  equal,  by  symmetry,  to  the  number  of 
antivortices,  this  result  predicts  that  the  number  of  vortices 
per  plaquette,  on  the  triangular  lattice,  in  the  infinite 
temperature  limit,  to  be  0.125.  Notice  that  this  number  is  in 
agreement  with  the  similar  result  derived  for  the  classical  XY 

model  on  this  lattice  in  section  A.  The  values  assumed  by 
2 

<p  >  and  <p'>  in  the  infinite  temperature  limit  are  0.125  and 
0.375  respectively.  Although  the  corresponding  limits  in  the 
s=°°  XY  model  have  not  been  calculated,  it  would  be  straight¬ 
forward  to  do  so.  What  can  be  said  in  this  regard  is  that  the 
number  of  vortex-antivortex  pairs  found  in  the  infinite 
temperature  limit,  on  both  lattices,  is  greater  than  that 
expected  by  assuming  the  probability  of  finding  a  pair  is  just 
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the  product  of  the  probability  of  locating  a  vortex  and  that 
of  locating  an  antivortex.  It  seems  that,  even  in  the  infinite 
temperature  limit,  the  chance  of  finding  an  antivortex  (or  a 
vortex  for  that  matter)  adjacent  to  a  vortex  is  enhanced 
relative  to  that  of  finding  one  of  these  structures  alone. 

This  is  actually  reasonable  when  it  is  noted  that,  on  the 
triangular  lattice  for  example,  a  vortex  configuration  on  a 
plaquette  involves  the  orientation  of  spins  at  three  lattice 
sites  while  the  addition  of  an  antivortex  to  the  adjacent 
plaquette  is  dependent  only  on  the  orientation  of  the  spin  on 
the  one  remaining  site  which  defines  the  domino  .  A  similar 
situation  occurs  on  the  square  lattice. 

In  the  ground  state  of  the  classical  (s=°°)  XY  model, 
which  consists  of  all  spins  aligned  in  some  arbitrary  direction 
in  the  x-y  plane,  there  are  no  vortices.  This  is  consistent 
with  the  view  that  the  vortices,  and  antivortices,  are 
topological  excitations  whose  numbers  grow  as  a  function  of  T. 
In  contrast  with  this  is  the  result  derived  above,  that  the 
number  of  vortices  and  antivortices  per  lattice  site  in  the 
ground  state  of  the  spin  h  XY  model  is  0.025+0.002  on  the 
square  lattice  and  0.104+0.008  on  the  triangular  lattice.  This 
is  to  be  considered  along  with  the  fact  that  the  ground  state 
of  this  model  shows  partial  long-range  order,  in  the  form  of  a 
non-zero  transverse  magnetization  per  site.  This  order  is  not 
total,  leaving  enough  local  disorder  to  allow  the  presence  of 
vortices  and  antivortices,  for  it  is  the  nature  of  these 
excitations  to  arise  as  local  magnetic  order  gives  way  to 
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disorder.  The  origin  of  both  this  disorder  and  the  non-zero 

vorticity  at  zero  temperature  is  of  course  quantum  mechanical, 

with  the  Hamiltonian  composed  of  the  non-commuting 

X  v 

operators  S^,  S  .  ,  so  that  the  ground  state  cannot  be  an 

eigenstate  of  either  operator.  It  is  interesting  to  note  that 

while  the  transverse  magnetizations  on  the  square  and  triangular 

lattices  differ  by  only  ^4%  (see  Marland  and  Betts,  {1.17}) 

the  vorticity  per  lattice  site  on  the  triangular  lattice  is 

almost  exactly  four  times  that  found  for  the  square  lattice. 

This  is  only  part  of  the  evidence  which  indicates  that 

understanding  of  the  combinatorics  of  vorticity  on  different 

lattices  is  far  from  complete. 

2 

Once  the  vorticity,  <v  >  and  <v' >  and  the  pair  vorticity, 

2 

<p  >  and  <p'>,  are  known  it  is  interesting  to  try  to  calculate 

the  number  of  free  vortices  present  in  the  ground  state  of  the 

spin  h.  XY  model.  By  "free  vortices "  what  is  meant  are  all  those 

vortices  (or  anti-vortices)  which  are  not  members  of  a  vortex- 

antivortex  pair  or  some  other  larger  cluster.  With  only  the 

total  vorticity  and  the  number  of  vortex-antivortex  pairs, 

2 

from  <p  >  and  <p'>,  available,  it  is  clear  that  the  effect  of  larger 
clusters  must  be  ignored  in  the  calculation  of  the  number  of 
free  vortices.  For  the  square  lattice,  where  the  density  of 
vortices  and  anti vortices  is  ^2.5%  this  contribution  is 
likely  negligible.  For  the  triangular  lattice  this  density  is 
^10%  and  the  effect  of  clusters  becomes  more  important.  This 
should  be  kept  in  mind  when  interpreting  the  results  presented 


below. 
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The  number  of  free  vortices  per  site  is  obtained  by 

associating  one  bound  vortex  with  every  pair  and  subtracting 

the  number  of  bound  vortices  from  the  total  number  of  vortices. 

The  latter  number  is  of  course  given  by  half  the  total  number 

2 

of  vortices  and  antivortices  ,  obtained  from  <v  >  and  <v*>. 

In  terms  of  operator  expectation  values  the  number  of  free 
vortices  per  site  on  the  square  lattice  is  then 

vfree  =  l<v2>  '  <p2>  '  (3'25) 


and  on  the  triangular  lattice  is 


v 


free 


> 


<P 


> 


(3.26) 


The  values  obtained  for  this  quantity  are  shown  in 

Table  3.1.  The  results  on  the  square  lattice  are,  with  one 

exception,  all  consistent  with  zero  in  the  sense  that  vfree<10 

2 

which  is  of  the  same  order  as  the  uncertainty  in  <v  >  and 
2 

<p  >  introduced  by  the  numerical  computations.  The  exception 

occurs  for  N=10  and  although  vfree  is  still  small  in  absolute 

2 

magnitude  it  is  of  the  same  order  as  <p  >^g.  ln  fact,  if  as 

2 

mentioned  above,  the  value  of  <p  >^g  were  doubled,  then  the 

number  of  free  vortices  per  site  for  this  cell  would,  with  a 

-4 

slight  extension  of  the  uncertainty  limit  to  v2xlo  ,  also  be 
consistent  with  zero.  The  results  for  the  triangular  lattice 
do  not  conform  as  well  with  expectation.  The  number  of  free 
vortices  per  site  is  both  relatively  large,  30%  of  <v*>/2,  and 
negative!  It  should  be  noted  however  that  in  absolute  terms 


71 


free  vortices  are  predicted  to  occupy  only  v-i%  of  the  lattice 

sites.  Possible  interpretations  of  and  explanations  for  this 

finding  will  be  speculated  upon  in  the  following  chapter. 

Finally,  it  is  interesting  to  compute  the  number  of 

free  vortices  per  site  for  the  infinite  lattice,  on  the  basis 

of  the  extrapolations  made  from  Figures  3.7  and  3.8.  For  the 

-4 

square  lattice,  the  calculation  yields  vfree=2x10  which,  in 

2 

view  of  the  large  uncertainty  of  0.001  in  <p  >^,13  consistent 
with  zero.  The  result  obtained  for  the  triangular  lattice  is 
vfree=-°-097±°.038  which  is  not  consistent  with  zero  free 
vortices  per  site,  even  with  the  inclusion  of  a  large 
uncertainty  in  <p,>co  ,  and  which  is  in  fact  not  physically 
possible. 

[  After  the  completion  of  the  main  body  of  this  work,  a 

new  formula  for  the  number  of  free  vortices  in  the  ground  state 
of  the  spin  h  XY  model  on  the  triangular  lattice  was  obtained 
{4.2}.  Let  p(v)  be  the  probability  of  finding  a  vortex  on  any 
given  plaquette  of  the  triangular  lattice  and  p(v;0)  the 
conditional  probability  that  an  adjacent  plaquette  contains 
neither  a  vortex  nor  an  antivortex.  Since  each  plaquette  on 
the  triangular  lattice  has  three  adjacent  neighbours,  the 
probability  of  finding  an  isolated  (or  free)  vortex  is  just 


p(v)  •  (p (v; 0) ) 


This  then  also  represents  the  number  of  free  vortices  present 
per  lattice  site.  On  the  triangular  lattice  p(v)=<v'>/2. 
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The  probability  p(v;0)  is  equal  to  l-p(v;a)  where  p(v;a)  is 
the  conditional  probability  that  a  plaquette  adjacent  to  one 
on  which  a  vortex  is  situated  is  itself  occupied  by  an  anti¬ 
vortex  (or  vortex) .  Neglecting  correlations  between  next- 
nearest  neighbour  plaquettes  (third  neighbour  and  higher  order 
site  correlations),  p(v;a)  is  just  the  probability  of  finding 
a  vortex-antivortex  pair  per  domino ,  <p'>/3,  divided  by  the 
probability  that  the  adjacent  site  is  occupied  by  an  antivortex, 
<  v 1  >/ 2  .  The  new  relation  which  gives  the  approximate  number 
of  free  vortices  is  then 


*  -  <v'  >  /  -i  2<p '  >  3 

vfree  2  3<v'> 


(3.27) 


Expanding  this  formula  to  first  order  in  <p’>/<v'>  gives 
equation  (3.26).  This  new  relation  is  an  improvement  over  the 
old  in  that  it  gives  a  direct  measure  of  the  number  of  free 
vortices,  although  it  still  neglects  the  effects  of  clusters 
of  three  or  more  vortices/antivortices  in  the  expression  for 

P  ( v ;  a )  . 

Using  the  values  of  <v'>  and  <p’>  in  Table  3.1b),  it  is 

k 

found  that  for  the  9-spin  cell  vfree=0 . 00691  while  for  the 

* 

13-spin  cell  vfree=0 • 00774 .  An  obvious  improvement  over  the 

k 

previous  values  for  the  number  of  free  vortices  is  that  vfree 
is  positive  for  both  values  of  N  for  which  data  is  available. 
These  numbers  are  also  smaller  in  absolute  value  than  the 
previous  estimates,  being  ^10%  of  the  respective  <v'>,  although 

k 

they  are  still  not  small  enough  to  be  consistent  with  vfree~0* 
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What  is  yet  to  be  calculated  is  the  expectation  value  of  the 
free  vortex  operator  which  is  introduced  in  Chapter  IV, 
section  B.  It  is  expected  that  such  a  quantity,  to  which 
(3.27)  is  only  an  approximation,  would  be  consistent  with  an 
absence  of  free  vortices  in  the  spin  %  XY  model  ground 
state.  ] 


' 


CHAPTER  IV 


SUMMARY  AND  DISCUSSION 

A.  Ground  State  Energy  and  Correlations 

The  main  numerical  result  of  Chapter  II  was  that  the 

ground  state  energy  of  the  spin  %  XY  model,  on  the  square  lattice, 

as  extrapolated  to  the  infinite  lattice  case  through  the 

results  of  calculations  on  cells  with  an  even  number  of  spins 

was  not  shown  to  be  equal  to  that  obtained  from  similar 

calculations  on  a  series  of  cells  containing  an  odd  number 

( 0 v0n ) 

of  spins.  The  corresponding  energies  are  Eq  /NJ  = 

-1.078  ±  0.002  and  E^^^/NJ  =  -1.07  ±  0.01.  The  uncertain- 

o 

ties  quoted  here  are  due  mainly  to  the  extrapolation  proce- 

-4 

dure,  with  computing  errors  accounting  for,  at  most,  5.0x10 
Certainly  before  any  conclusion  can  be  reached,  the  ground 
state  energy  for  N=17,  which  is  the  last  odd  number  of  sites 
per  cell  for  which  calculations  are  now  feasible,  must  be 
obtained.  The  results  of  the  even  N  calculations  indicate 
that  the  N=17  data  point  should  not  lie  more  than  ±0.004  units 
from  the  line  already  obtained  from  the  two  odd  N  data  points. 

Such  a  deviation  is  more  than  enough  however  to  bring  the  two 
extrapolated  ground  state  energies  into  agreement.  Although 
such  a  course  may  be  premature  it  is  interesting  to  speculate 
on  any  disagreement  between  the  infinite  lattice  energies 
E^-Vnu  and  E^°^^^/NJ  remaining  after  the  N=17  result 
has  been  calculated. 

The  source  of  any  discrepancy  may  be  due  to  the  fact 
that  the  periodic  boundary  conditions  force  a  tiling  of  the 
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infinite  square  lattice  with  domains  of  m  -h.  If  just  two 

z 

such  cells  of  dimension,  e.g.,  9  are  considered,  the  result 

is  a  configuration  of  18  spins  with  a  total  m of  1,  which 

is  taken  to  contribute  to  the  ground  state  of  the  infinite 

lattice.  The  argument  of  Mattis,  as  presented  in  Chapter  II  , 

would  imply  that  for  this  18  spin  block  the  ground  state 

has  m  =0.  An  extra  "domain  wall"  energy  has  been  added  that 
z 

does  not  disappear  in  the  infinite  lattice  limit,  thus  raising 

the  ground  state  energy  calculated  from  cells  with  an  odd 

number  of  spins  relative  to  that  obtained  from  even  N  cells. 

If  this  is  the  case,  better  agreement  may  be  found  by  taking 

"anti-periodic"  boundary  conditions  for  odd  N  cells,  in  which 

spins  on  sites  in  the  cell  boundary  are  linked,  not  to  other 

spins  on  boundary  sites,  but  to  the  opposite  spins  on  the  same 

boundary  sites.  Such  a  procedure  would  change  the  symmetry 

properties  of  the  cell,  but  it  would  have  the  effect  of 

patching  together  cells  with  opposite  net  m  ,  to  build  up 

an  infinite  lattice  which  at  each  stage  has  the  correct  m  , 

z 

according  to  Mattis'  arguments,  of  the  ground  state. 

The  lack  of  a  reliable  extrapolation  for  the  two  spin 
correlations  on  the  square  lattice  presented  in  Chapter  II 
has  already  been  mentioned.  To  eliminate  the  degeneracy  of 
the  second  and  third-nearest  neighbour  correlations  for  cells 
with  N£l6  which  seems  to  be  creating  the  difficulties,  one 
might  believe  that  performing  calculations  on  cells  with  less 
symmetry  than  the  square  cells  considered  here  would  be 
effective.  Such  calculations  were  however  carried  out  early 
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in  the  development  of  the  finite  cell  method.  The  results 
showed  no  pattern  at  all.  This  would  seem  to  indicate  that  a 
family  of  cells  which  share  the  same  symmetry  properties  is  an 
important  factor  in  achieving  linearity  with  this  method.  The 
four  spin  correlations  presented  in  Chapter  II  are  even  more 
erratic  as  functions  of  N  and  no  reliable  extrapolations  can 
be  obtained  for  these.  The  problem  here  seems  to  be  related 
to  cell  size  and  until  such  time  as  calculations  with  cells 
of  N>18  are  feasible  other  avenues  must  be  explored.  Once 
the  question  of  the  ground  state  energy  for  cells  with  an  odd 
number  of  spins  is  settled,  correlations  may  be  calculated  for 
these  cells.  The  extra  data  points,  while  not  as  useful  as 
results  for  N>18,  would  likely  give  an  indication  as  to  the 
direction  in  which  extrapolation  is  most  reliable.  As  for 
the  correlations  on  the  triangular  lattice  the  above  remarks 
may  certainly  apply,  but  the  immediate  concern  is  to  generate 
the  four  spin  correlations  for  N=19  and  see  whether  or  not 
the  results  form  a  straight  line,  before  apologizing  for  the 
lack  of  linearity  and  a  reliable  extrapolation. 

To  close  this  section  a  comment  on  the  applicability 
of  the  method  described  in  Chapter  II  to  other  models  besides 
those  already  considered  { 1 . 15 }- { 1 . 17 }  is  in  order.  An  example 
of  such  a  model  is  the  Ising  model  in  a  transverse  field, 
described  by  the  Hamiltonian: 


H  =  -  J 


<i  ,J> 


szsz  -  r 

i  D 


£  S * 

l 


l 


(4.1) 


. 


I  ‘ 
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This  model  is  particularily  interesting  since  Suzuki  has  proven 
{4.1}  the  equivalence  of  the  ground  state  of  (4.1)  with  that 
of  the  three-dimensional  Ising  model,  at  all  finite  tempera¬ 
tures.  Unfortunately,  this  model  is  not  amenable  to  analysis 
by  the  finite  cell  method  since  there  appears  to  be  no  simple 

many -body  operator  such  as  Mz  =  E  S.,  which  commutes  with  H 

i 

and  thus  provides  a  good  quantum  number  to  label  the  energy 

states,  in  particular  the  ground  state.  The  entire  space  of 

N 

states,  of  dimension  2  for  a  cell  of  N  spin  V s ,  would  seem 
to  contribute  to  the  ground  state  and  for  N^IO,  even  taking 
the  symmetry  of  the  cell  into  account,  computations  are  at 
present  unfeasible.  The  XY  model  in  a  longitudinal  field 
does  not  suffer  from  this  defect,  since  M  still  commutes 

Zj 

with  the  Hamiltonian. 

H  =  -2 J  E  (SXSX  +  sYsY)  -  r  E  s2  ,  (4.2) 

<i,j>  13  1  :  i  1 

and  although  m  is  not  necessarily  zero  in  the  ground  state 

z 

for  a  certain  range  of  the  parameter  T/J,  the  calculations 
involved  are  no  more  complex  than  those  already  carried  out 
for  H  .  Such  calculations  remain  to  be  performed. 

X  Y 

B .  Vortices  and  Free  Vortices 

The  results  of  the  calculation  of  the  total  number  of 

vortices  and  antivortices  per  lattice  site  in  the  ground 

2 

state  of  the  spin  h  XY  model  are  <v  >^=0.025  ±  0.002  on  the 
square  lattice  and  <v'>oo=0.104  ±  0.008  on  the  triangular 
lattice.  This  is  unlike  the  situation  in  the  ground  state 
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of  the  classical  (s=°°)  XY  model,  which  contains  no  vortices 

or  antivortices.  The  interpretation  of  these  findings  in 

terms  of  partial  long-range  order  and  local  disorder  has 

already  been  made  in  the  previous  chapter.  A  comment  also 

made  previously  but  which  bears  repeating  deals  with  the 

relative  magnitudes  of  the  total  vorticity  on  the  square  and 

triangular  lattices.  Although  it  is  certainly  not  a  definite 

indication  of  any  relation,  the  fact  that  the  triangular 

lattice  result  is  almost  exactly  four  times  the  square  lattice 

result  suggests  a  geometric  interpretation  may  be  possible 

for  this  difference,  especially  in  view  of  the  fact  that  the 

transverse  magnetization,  which  signals  long-range  order,  is 

so  similar  on  both  lattices.  Another  integer  factor  arises 

when  the  expectation  value  of  the  pair  operator  for  the  N=10 

cell  on  the  square  lattice  is  considered.  The  individual 

2 

two  and  four  spin  correlations  which  contribute  to  <p  >^q 

do  not  seem  to  lie  outside  of  the  expected  range,  but  an 

examination  of  Figure  3.8  as  well  as  consideration  of  the  free 

vortices  per  site  for  the  N=10  cell  indicates  that  a  value  of 
2 

<p  >^q  twice  that  obtained  would  be  in  better  agreement  with 

the  other  results  on  the  square  lattice.  The  only  place  for 

such  a  factor  to  enter  is  through  the  geometry  of  the  cell. 

There  is  a  clue  here  as  well,  for  the  N=10  cell  is  the  only 

one  considered  on  the  square  lattice  which  was  not  reflection 

invariant.  The  exact  effect  this  might  have  on  the  calculation 
2 

of  <p  >-^Q  is  unknown,  but  it  does  suggest  that  cell  geometry 
plays  a  larger  role  than  assigned  in  the  determination  of 
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vorticity  and  pair  vorticity. 

In  the  classical  (s=°°)  XY  model,  in  two  dimensions  of 
course,  the  introduction  of  vortices  and  antivortices  as 
excitations  of  the  system  has  allowed  an  intuitive  picture 
of  the  phase  transition  which  has  been  predicted  for  this 
model  to  be  developed.  Vortex-antivortex  pairs  are  created 
as  temperature  rises,  until,  at  the  critical  temperature, 
these  pairs  unbind  (the  vortex-antivortex  "dipoles"  are 
"ionized")  and  the  long  range  order  which  characterizes  the 
low  temperature  phase  is  lost.  A  phase  transition  has  also 
been  predicted  in  the  spin  %  XY  model  and  so  an  interesting 
quantity  to  calculate  in  the  ground  state  of  this  model  is 
the  number  of  free  vortices,  those  which  are  not  found 
adjacent  to  an  anti vortex.  Although  the  number  of  vortices 
in  the  ground  state  of  the  spin  h  XY  model  is  non-zero,  the 
number  of  free  vortices  is  expected  to  be  zero. 

Results  of  the  calculation  of  the  number  of  free 
vortices  in  the  ground  state  of  the  spin  h  XY  model  on  the 
square  lattice  confirm,  with  the  one  well-noted  exception 
discussed  above,  the  predicted  value  of  zero.  This  result 
holds  in  the  infinite  lattice  limit  as  well.  The  number  of 
free  vortices  on  the  triangular  lattice,  as  derived  above,  is 
not  consistent  with  zero  and  is  in  fact  negative,  a  result 
which  holds  in  the  infinite  lattice  limit  as  well.  The  origin 
of  this  result  seems  to  be  the  larger  density  of  vortices  and 
antivortices  which  occurs  in  the  ground  state  on  the  triangu- 
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lar  lattice  relative  to  the  square  lattice.  The  numbers  are 
10.4%  relative  to  2.5%.  When  this  density  increases,  the 
number  of  pairs  which  are  free  decreases  as  larger  clusters 
of  three  and  four  vortices/antivortices  are  formed.  When 
this  occurs  it  is  no  longer  true  that  one  vortex  is  associated 
with  every  vortex-antivortex  pair  as  for  example  a  cluster  of 
two  vortices  and  two  antivortices  on  the  triangular  lattice 
will  contain  3  such  pairs.  The  number  of  bound  vortices  will 
be  overestimated  and  the  number  of  free  vortices  obtained 
from  the  relation  (3.23)  may  even  be  negative.  As  an  example 
consider  the  anti-ferromagnetic  ground  state  of  the  classical 
(s=o&)  XY  model  on  the  triangular  lattice  (see  Appendix  III)  . 

In  this  case  the  number  of  pairs  per  domino  is  1,  since  the 
state  is  "packed"  with  vortices  and  antivortices  and  every 
domino  contains  a  pair.  The  number  of  pairs  per  site  is  3. 

The  number  of  vortices  per  plaquette  is  i,  every  site  being 
occupied  by  either  a  vortex  or  antivortex,  which  occur  in 
equal  numbers,  so  that  the  number  of  vortices  per  site  is  1. 
Using  (3.23)  to  calculate  the  number  of  free  vortices,  one 
finds  this  number  to  be  -21  The  reason  for  this  should  now 
be  apparent.  Just  as  there  really  are  no  free  vortices  in 
this  state,  there  are  no  free  pairs  either.  Every  pair  is  in 
turn  a  member  of  some  higher  order  cluster  which  in  turn  is 
part  of  an  even  larger  cluster,  and  so  on  ad  infinitum.  There 
are  no  sites  in  the  lattice  for  this  state  which  are  not 
occupied  by  a  vortex  or  antivortex.  The  simple  formula  (3.23) 
is  not  applicable  then  in  situations  where  the  total  number  of 
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pairs  and  the  total  number  of  free  pairs  are  not  the  same. 

It  would  then  seem  reasonable  to  conjecture  that  the 
fact  that  equation  (3.23)  predicts  a  negative  number  of  free 
vortices  on  the  triangular  lattice  is  an  indication  that  the 
number  of  three  and  especially  four-member  clusters  of  vortices 
and  antivortices  is  non-zero  in  the  ground  state  of  the  spin 
h  XY  model  on  the  triangular  lattice.  This  result  also 
suggests  that  the  number  of  free  vortices  is  indeed  zero  in 
the  ground  state  on  this  lattice. 

Verification  of  this  can  be  sought  in  the  calculation 
of  the  expectation  value  of  the  free  vortex  operator  which 
has  recently  been  calculated  {4.2}  for  the  triangular  lattice. 
The  eigenstates  of  this  operator  are  taken  to  be  those  in 
which  a  vortex  situated  on  a  central  plaquette  is  surrounded 
on  the  three  adjacent  plaquettes  by  configurations  exhibiting 
no  vorticity.  A  total  of  six  sites  are  involved  and  the 
expectation  value  of  this  operator  involves  all  three  two 
spin  correlations  (first,  second  and  third-nearest  neighbour) , 
all  inequivalent  four  spin  correlations  and  several  six  spin 
correlations  .  Pressures  of  time  have  not  allowed  the 
evaluation  of  this  operator,  although  it  is  straightforward 
to  accomplish.  It  should  be  noted  however  that  because 
correlations  between  spins  up  to  third-nearest  neighbour  in 
distance  must  be  evaluated,  calculations  would  be  required  on 
the  19- spin  cell  in  order  to  achieve  an  accurate  extrapolation 
to  the  infinite  lattice.  Since  in  fact  the  9  and  13-spin  cells 
do  not  have  third  nearest  neighbours,  a  simpler  free  vortex 


82 


operator  can  be  derived  which  is  applicable  to  these  cells. 

In  principle  the  same  operator  could  be  determined  for  the 
square  lattice.  A  number  of  reasons  suggest  this  would  be 
of  little  use.  First,  the  number  of  free  vortices  in  the 
ground  state  of  the  spin  h  XY  model  on  the  square  lattice 
has  been  estimated  to  be  zero  using  the  approximate  relation 
(3.22).  Secondly  the  construction  of  such  an  operator  would 
involve  consideration  of  configurations  over  twelve  sites 
and  although  possible,  it  would  be  likely  that  computer 
assistance  would  be  necessary  to  complete  the  derivation. 
Finally,  such  an  operator  would  involve  correlations  between 
spins  separated  by  a  distance  of  up  to  /lO  6,  and  it  is  likely 
that  the  numerical  results  for  the  16  and  18  spin  cells, 
the  only  two  for  which  this  operator  expectation  value  exists, 
would  not  be  enough  to  allow  for  an  accurate  extrapolation 
to  the  infinite  lattice. 

C .  Experiments  on  Two-Dimensional  Systems 

The  discussions  to  this  point  have  all  been  of  a 
theoretical  character,  with  little  mention  of  experiment. 

The  Kosterlitz-Thouless  theory  has  been  taken  as  the 
description  of  the  phase  transition  in  the  planar  model  and 
competing  theories  have  been  neglected.  The  reason  this 
approach  can  be  taken  is  that  experimental  evidence  now  favors 
the  Kosterlitz-Thouless  picture  over  those  of  other  investi¬ 
gators.  To  discuss  these  experiments  however,  one  final 
theoretical  result  is  needed.  Applying  the  Kosterlitz- 
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Thouless  theory  to  neutral  superfluids.  Nelson  and  Kosterlitz 

{4.3}  discovered  a  simple  relationship  between  superfluid 

density,  p  ,  and  the  critical  index  r\ (T) . 
s 
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'2  7  7 

h  n(T) 


(4.3) 


where  m  is  the  mass  of  a  Helium  atom,  h  is  Planck's  constant 
and  the  index  r\  (T)  is  defined  from  the  spin-spin  correlation 
function  by 
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(4.4) 


for  T<T  .  The  relation  (4.3)  holds  in  the  limit  as  T->T 
c  c 

from  below.  Thus  a  measurement  of  the  jump  in  the  superfluid 

density  p  at  the  critical  temperature  can  be  compared  with 
s 

the  value  predicted  from  (4.3)  for  this  quantity: 

2 

p  (T)  2nm  k 

lim_  {-| - }  =  .  (4.5) 

T+T  h  n (T  ) 

c  c 

The  Kosterlitz-Thouless  theory  (1.14),  {4.4}  predicts  the 

universal  result  ri(T  )=^  which  leads  from  (4.5)  to 

P  c  ( 1 )  _  q  9 

lim_  {— - - }  =  3.491x10  g/cm  -K 

T+T  1 
c 

Several  experiments  have  now  been  carried  out  which 

measure  the  superfluid  density  at  the  transition.  Rudnick 

and  collaborators,  beginning  in  1969,  had  performed  measure- 

4 

ments  of  third  sound  propagation  in  thin  He  films  {4.5} 
and  a  reanalysis  of  this  earlier  data  yields  {4.6} 
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(3.30  ±  0 . 21) xlO  9  g/cm2-K 


Bishop  and  Reppy  {4.7}  have  measured  the  superfluid  transi- 


4 

tion  in  two-dimensional  He  films  by  "shaking"  the  fluid, 


using  the  Andronikash vili  torsional-oscillator  technique. 

A  helium  film  is  adsorbed  on  a  Mylar  film  substrate.  A 
strip  is  wound  up  into  a  cylinder  and  oscillated  at  the 
resonant  frequency.  When  a  helium  film  of  the  desired 
thickness,  usually  one  or  two  monolayers,  has  been  deposited, 
measurements  of  the  period  and  Q  factor  of  the  oscillator 
are  made  to  determine  the  temperature  dependence  of  the 
superfluid  mass  and  dissipation.  The  value  obtained  for 


p  IT  ) /T  is 
s  c  c 


These  experiments  give  strong  support  to  the  Kosterlitz- 
Thouless  version  of  the  planar  critical  behaviour  and  also 
appear  to  eliminate  several  competing  theories  which  predict 
other  values  for  n (T  ) . 


c 


Other  systems  of  experimental  interest  include  {4.8}, 
the  solid  formed  by  electrons  on  the  surface  of  liquid  helium, 
at  sufficiently  low  temperatures  and  high  electron  density, 
various  systems  liquid  crystals,  lipid  layers  floating  on  a 
water  surface  and  layers  adsorbed  on  a  crystalline  substrate. 
All  of  these  systems  can  be  discussed  in  terms  of  the  theory 
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of  two-dimensional  melting  {l.ll},  with  various  modifications, 

but  at  the  moment  no  firm  connection  between  theory  and 

experiment  has  been  made.  In  fact  the  electron  solid  has 

only  recently  been  observed  {4.9},  using  a  radio-f requency 

resonance  technique.  Another  recent  discovery  is  the 

construction  of  an  physical  system  which  is  apparently  a 

two-dimensional  spin  h  XY  antiferromagnet .  The  system 

consists  of  crystals  of  (Co (p-methyl-pyridine  N-oxide)^) 

(CIO^^  (4.10),  with  the  Co  ion  in  a  high  ligand  field 

possessing  spin  At  finite  (T^T  and  above)  temperatures, 

c 

the  spin  %  XY  model  is  also  applicable  to  a  variety  of 

magnetic  insulators  {4.11}  ,  although  experimental  systems 

which  remain  well  described  by  the  XY  model  at  very  low 

(T<<T  )  temperatures  have  not  yet  been  discovered.  It 
c 

remains  very  much  an  experimental  challenge  to  provide 
results  which  will  complement  theoretical  findings  on  the 
nature  of  the  phase  transition  in  the  two-dimensional  spin  h 


XY  model. 
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APPENDIX  I 


TWO  SPIN  CORRELATIONS  IN  THE  GROUND  STATE  OF  THE 
XY  MODEL  ON  THE  SQUARE  LATTICE  FINITE 

CELL  RESULTS 

In  Chapter  II  only  those  two  spin  correlations  which 
were  needed  for  the  calculation  of  the  vorticity  and  pair 
vorticity  operators  were  presented.  A  large  number  of  such 
correlations  have  been  calculated  however  and  Table  AI . 1  is  a 
compilation  of  both  transverse  and  longitudinal  correlations 
for  finite  cells  on  the  square  lattice.  The  correlations  for 
even  N  cells  were  calculated  by  Dr.  J.Oitmaa  ,  as  part  of  the 
research  described  in  {1.15}  while  those  for  odd  N  cells  were 
calculated  as  part  of  the  research  described  in  this  thesis. 

The  correlations  are  defined  as 

C_  (r/S )  =  <SXSX>  =  t<S+S-  +  sV>  ,  (AI.l) 

T  or  4or  ro 

and  CT  (r/6)  =  <SZSZ>  ,  (AI.2) 

L  or 

where  S  is  a  spin  h  operator  and  6  is  the  lattice  spacing.  The 

values  of  r/6  listed  in  the  table  are  those  for  which  these 

correlations  are  independent.  For  example,  for  N=10,  each  site 

has  two  distinct  first  and  second-nearest  neighbours,  while 

there  are  no  independent  third-nearest  neighbour  pairs  and  only 

16  fourth-nearest  neighbour  pairs.  For  N=8,  there  are  only 

first,  second  and  third-nearest  neighbours,  16  pairs  of  each. 

-5 

Five  figures  are  included  and  an  uncertainty  of  ±  1.0x10  due 
to  machine  errors  may  be  assumed. 
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Table  AI . 1  -  Transverse  and  longitudinal  correlations  in  the 

ground  state  of  the  spin  h  XY  model  on  finite 
cells  of  the  square  lattice. 


N 

r/6 

CL(r/S) 

CT(r/«) 

4 

1 

-0.12500 

0.17678 

v/2 

0 

0.12500 

5 

1 

-0.05000 

0.15000 

8 

1 

-0.05398 

0.14657 

/2 

-0.01136 

0.13078 

2 

-0.01136 

0.13078 

9 

1 

-0.04652 

0.14371 

vT 

-0.00904 

0.12493 

10 

1 

-0.05029 

0.14402 

/2 

-0.00958 

0.12686 

/T(  16) 

-0.01050 

0.12192 

13 

1 

-0.04440 

0.14070 

ST 

-0.00707 

0.12393 

2 

-0.00622 

0.12105 

16 

1 

-0.04558 

0.19062 

/2 

-0.00663 

0.12171 

2 

-0.00663 

0.12171 

/5 

-0.00585 

0.11696 

/ 8 

-0.00444 

0.11445 
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Table  AI . 1  -  Transverse  and  longitudinal  correlation  in  the 

ground  state  of  the  spin  %  XY  model  on  finite 
cells  of  the  square  lattice.  (Continued) 


N 


r/6 


CL(r/6)  CT(r/6) 


1 

-0.04445 

0.13988 

/T 

-0.00629 

0.12196 

2 

-0.00513 

0.11819 

/r 

-0.00542 

0.11744 

3 

-0.00486 

0.11590 

APPENDIX  II 


GXYSQDATA  and  GXY4C0R 

In  this  brief  appendix  the  content  of  the  files 
GXYSQDATA  and  GXY4C0R,  which  have  been  placed  in  the 
Statistical  Physics  Program  and  Data  Library,  will  be 
discussed.  Only  a  description  of  the  contents  of  these 
files  is  presented.  Details  can  be  obtained  from  the  files 
themselves,  which  have  been  documented  for  this  purpose. 

The  file  GXYSQDATA  contains  all  the  data  which  has 
been  generated  by  finite  cell  calculations  on  the  square 
lattice.  For  N=8,9,10  and  13  this  includes  the  elements  of 
the  space  group  S  ,  the  basis  states  and  their  symmetry 
numbers,  the  matrix  elements  of  -H  /J  in  this  basis  and  the 

A.  x 

eigenvalues  and  eigenvectors  of  this  matrix.  Because  of  the 
large  size  of  the  Hamiltonian  matrix  for  N=16  and  18,  only 
the  greatest  eigenvalue  and  the  corresponding  eigenvector 
are  recorded  for  these  cells.  In  addition  the  limitations  of 
storage  space  prevented  the  Hamiltonian  matrix  elements  for 
the  N=18  cell  from  being  written  into  a  permanent  file.  These 
therefore  are  not  present  in  GXYSQDATA.  A  list  of  nearest 
neighbour  pairs,  for  each  N,  is  also  present  in  this  file. 
Together  with  a  knowledge  of  the  cell  geometry,  this  data 
allows  the  numbering  of  each  cell  which  was  used  to  be 
reproduced.  This  knowledge  is  important  for  the  interpreta¬ 
tion  of  the  elements  of  and  the  basis  states  in  geometric 
terms.  The  line  format  of  each  piece  of  data  is  recorded  in 
the  file. 
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The  program  file  GXY4C0R  contains  the  program  written 
to  evaluate  four  spin  correlations  of  the  form 


<^N)  |  (il,jk)  +  (ik, jl)  -  (ij,kl)  |«^N)> 


(All. 1) 


where 


(i j  , kl)  =  s!s+S,  S-. 

1  j  k  1 


(All .2) 


The  dimension  and  format  statements  included  are  appropriate 
for  calculations  on  16-spin  cell  data.  As  input,  the  program 
requires  the  elements  of  S  ,  the  basis  states  and  their 
symmetry  numbers  and  the  coefficients  of  the  ground  state 
eigenvector  in  this  basis.  This  data,  for  all  N,  is  present 
in  GXYSQDATA.  For  any  N£l3,  the  program  is  designed  to  calcu¬ 
late  up  to  five  correlations  of  the  form  (AII.l).  The  numbers 
i,j,k  and  1  are  specified  directly  in  the  program  and  need  not 
be  included  in  the  input  file.  For  N=16  limitations  of  central 
core  space  which  is  available  to  an  individual  program 
prevents  more  than  two  such  correlations  from  being  evaluated, 
during  one  run  of  the  program.  For  N=18  and  19  (the  latter 
being  the  dimension  of  a  cell  on  the  triangular  lattice)  only 
one  correlation  of  the  form  (AII.l)  can  be  computed  during 
one  run.  An  interesting  feature  of  the  program  is  that,  just 
as  it  allows  the  direct  calculation  of  (AII.l)  which  involves 
a  sum  (and  difference)  of  correlations  of  the  form 


any  ground  state  correlation  of  the  type 
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<$^N)  |a  (ij,kl)  +  a2(i,j',k,l')  +  a  (i"j",k"l")  + 


$<N>> 

o 


where  a^,a2,a^  are  any  numbers  and  the  sets  {i,j,k,l}, 

{ i ' ,k' ,k' , 1* }  and  {i" , j " ,k" , 1" } , . . .  are  all  distinct,  can  be 
calculated  with  small  modification  to  GXY4C0R.  In  the  case 
of,  e.g.  N=19 ,  this  feature  will  allow  the  direct  calculation 
of  the  sum  of  the  two  four  spin  correlations  involved  in  <p'> 
in  approximately  the  same  amount  of  time  as  the  computation 
of  one  such  correlation  would  require.  Because  of  the  large 
amounts  of  time  involved  in  these  calculations  such  a  saving 
may  outweigh  any  advantage  in  calculating  both  correlations 
involved  separately. 


APPENDIX  III 


ON  THE  GROUND  STATE  OF  THE  ANTIFERROMAGNETIC 

PLANE  ROTATOR  MODEL 

By  analogy  with  the  ferromagnetic  plane  rotator  model 
introduced  in  Chapter  I,  the  antiferromagnetic  plane  rotator 
is  defined  by  the  Hamiltonian 

H  =  J  E  S.*S.  =  J  E  cos  (<J> .  —  4>  . )  ,  (AIII.l) 

.  .  11  ..  l  i 

<i,D>  3  <i / 3  >  J 

where  the  coupling  constant  J  is  positive  and  <J> .  is  the  angle 
made  by  the  spin  vector  with  some  reference  axis.  The 

■V  .  #  (  t  , 

vector  is  assumed  to  have  unit  magnitude.  The  minimum 
energy  configuration  of  spins  which  interact  according  to 
this  Hamiltonian  is  such  that  neighbouring  spins  and 
are  anti-parallel,  or  as  close  to  anti-parallel  as  possible. 
Now  both  the  square  and  hexagonal  lattices  can  be  divided 
into  two  sub-lattices,  A  and  B,  such  that  all  the  nearest 
neighbours  of  a  particular  spin  on  one  sub-lattice  belong  to 
the  other  sub-lattice.  The  antiferromagnetic  ground  state  is 
then  clearly  that  state  in  which  all  the  spins  on  sub-lattice 
A  are  oriented  in  one  direction  and  all  the  spins  on  sub¬ 
lattice  B  are  oriented  in  the  opposite  direction.  Each 
nearest  neighbour  interaction  is  then  -1  and  the 

resulting  energies  per  site,  -2J  on  the  square  lattice  and 
-3/2J  on  the  hexagonal  lattice,  are  minimal.  Note  that  there 
is  a  degeneracy  in  this  type  of  state,  in  the  sense  that  the 
angle  of  the  spins  on  sub-lattice  A  with  respect  to  a  fixed 
reference  axis  can  assume  any  value  in  the  range  0  to  2tt  rad. 
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This  degeneracy  is  of  course  shared  with  the  ferromagnetic 
version  of  (AIII.l). 

The  same  procedure  cannot  be  applied  to  the  triangular 
lattice  which  can  in  fact  be  divided  into  a  minimum  of  three 
geometrically  similar  sub-lattices.  This  lattice  however 
may  be  regarded  as  constructed  out  of  unit  triangular  cells 
and  the  construction  of  the  ground  state  begins  with  the 
determination  of  the  minimum  energy  configuration  of  one 
such  cell.  With  the  cell  vertices  numbered  1,  2  and  3,  and 
taking  the  spin  on  site  1  to  define  a  reference  axis,  the 
cell  Hamiltonian  can  be  written  as 

H  n  =  J  ( cos  0  0  +  cos  0  _  +  cos(6„-6-j))  ,  (AIII.2) 

cell  2  3  2  3 

where  the  angles  ©2  and  ©^  are  measured  with  respect  to  the 
reference  axis.  Minimizing  this  expression  as  a  function  of 
©2  and  0^  leads  to  the  following  two  distinct  solutions  for 
02  and  0^: 

=  “2  it/ 3  ,  e*1’  =  +2tt/3  ,  (AIII.3) 

and  0*2)  =  +2ir/3  ,  e.*2^  -  2ir/3  .  (AIII.4) 

With  a  positive  angle  corresponding  to  counter-clockwise 
rotation,  and  with  a  numbering  of  vertices  as  in  Figure  3.4, 
it  is  easy  to  see  that  the  solution  B^^ /  B^"^  represents 
the  vortex  configuration  of  Figure  3.4  (a) 1.  While  the 
solution  6^  /  03  corresponds  to  the  antivortex  configura¬ 

tion  of  Figure  3.4(b)l.  Again  there  is  a  degeneracy  implicit 
in  the  two  configurations  represented  by  these  solutions,  in 
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the  sense  that  the  orientation  of  the  spin  on  site  1  is 
arbitrary.  Given  that  the  spin  one  one  site  of  a  unit  tri¬ 
angular  cell  has  an  arbitrary,  fixed  direction,  two  config¬ 
urations  that  have  the  same,  minimum,  energy  can  be  construc¬ 
ted.  One  is  a  vortex  and  the  other  is  an  antivortex.  These 
are  generated  by  first  numbering  the  vertices  of  the  cell  so 
that  the  positive  direction  coincides  with  123  (as  in  Figure 
3.4),  where  site  1  is  the  location  of  the  fixed  spin.  A 
vortex  is  constructed  on  this  cell  by  placing  the  spin  on  site 
2  at  an  angle  of  -2tt/3  rad.  with  respect  to  the  fixed  spin 
on  site  1,  and  the  spin  on  site  3  at  the  same  angle  of  -2tt/3 
rad.  but  now  with  respect  to  the  spin  on  site  2.  The  original 
spin  on  site  1  is  "automatically"  oriented  at  an  angle  of 
-2tt/3  rad.  with  respect  to  the  spin  on  site  3.  An  antivor¬ 
tex  is  constructed  by  following  the  same  procedure  with  a 
displacement  angle  of  +2u/3  radians. 

It  is  now  possible  to  construct  the  ground  state  of 
the  antiferromagnetic  plane  rotator  on  the  triangular  lattice 
and  understand  why  it  is  the  minimum  energy  state.  First 
divide  this  lattice  into  three  sub-lattices,  A,  B  and  C  as 
shown  in  Figure  AIII.l.  Place  all  the  spins  on  sub-lattice  A 
parallel  to  each  other  and  oriented  in  some  arbitrary  direc¬ 
tion.  Then  place  all  the  spins  on  sub-lattice  B  in  a 
direction  rotated  by  -2tt/3  rad.  with  respect  to  the  orienta¬ 
tion  of  the  spins  on  the  A  sub-lattice.  The  spins  on  the  C 
sub-lattice  are  then  oriented  at  an  angle  of  -2tt/3  rad.  with 
respect  to  the  spins  on  the  B  sub-lattice.  Now  each  triangular 

cell  on  the  lattice  belongs  to  one  of  six  classes,  with 
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respect  to  the  orientation  of  the  site  spins  as  each  cell 
is  traversed  in  a  positive  direction;  these  are  ABC,  BCA, 

CAB,  ACB ,  CBA  and  BAC .  The  first  three  and  last  three  are 
easily  seen  to  correspond  to  degenerate  energy  states,  since 
the  energy  of  a  unit  triangular  cell  depends  only  upon  the 
relative  orientation  of  the  spins  about  this  cell.  The  form 
ABC  corresponds  to  a  vortex  by  construction  so  the  form  ACB 
must  correspond  to  an  antivortex,  as  can  be  seen  from 
Figures  3.4(a)l  and  (b)  1.  Now  it  was  shown  above  that  the 
vortex  and  antivortex  configurations  are  degenerate  in  energy 
and  that  this  energy  was  the  minimum  possible  for  a  unit 
triangular  cell.  Because  each  and  every  such  cell  on  the 
lattice  of  Figure  AIII.l  is  in  the  minimum  energy  configura¬ 
tion,  this  entire  lattice  is  in  a  minimum  energy  configuration. 

The  ground  state  so  derived  is  illustrated  in  Figure 
AIII.2  for  a  particular  choice  of  the  direction  of  the  spins 
on  sub-lattice  A.  It  is  interesting  to  note  that  this  state 
consists  of  vortices  and  antivortices  packed  into  the  lattice 
so  that  each  vortex  (antivortex)  is  adjacent  to  three  anti¬ 
vortices  (vortices).  No  plaquette,  in  the  terminology  of 
Chapter  III,  has  zero  vorticity.  Another  interesting  feature 
of  this  state,  mentioned  in  Chapter  III,  is  that  if  it  is 
regarded  as  constructed  out  of  "upright"  triangular  cells 
with  vortex  configurations  of  spins,  antivortex  configurations 
will  appear  in  the  "inverted"  triangular  cells  which  appear 
as  the  "upright"  cells  are  joined  together  to  form  the 
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Figure  AIII.l  -  Triangular  lattice  partitioned  into  three 

geometrically  similar  sub-lattices. 


Figure  AIII.2  -  Antiferromagnetic  ground  state  of  the  plane 

rotator  model  on  the  triangular  lattice 
(vEvortex,  a Eanti vortex) . 
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infinite  lattice.  It  thus  appears  that  the  triangular 
lattice  cannot  be  tiled  with  vortex  configurations  alone 
without  generating  antivortices  as  well,  although  the 
geometry  of  the  vortex  and  antivortex  cells  is  slightly 
different,  as  noted  in  Chapter  III. 


